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Abstract: 

We propose to resum exactly any number of one-loop vacuum polarization insertions 
into the scale of the coupling of lowest order radiative corrections. This makes maximal 
use of the information contained in one-loop perturbative corrections combined with 
the one-loop running of the effective coupling and provides a natural extension of the 
familiar BLM scale-fixing prescription to all orders in the perturbation theory. It is 
suggested that the remaining radiative corrections should be reduced after resummation. 
In this paper we implement this resummation by a dispersion technique and indicate 
a possible generalization to incorporate two-loop evolution. We investigate in some 
detail higher order perturbative corrections to the r decay width and the pole mass of 
a heavy quark. We find that these corrections tend to reduce a s (m T ) determined from 
r decays by approximately 10% and increase the difference between the bottom pole 
and MS-renormalized mass by 30%. 
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1 Introduction 



During the past decade QCD has turned from qualitative descriptions to quantita- 
tive predictions for the strong interactions. The influx of more accurate experimental 
data represents a constant stimulus to improve theoretical predictions. While an un- 
derstanding of confinement remains as evasive as ever, these predictions rely on the 
validity of nonperturbative factorization in processes governed by a large momentum 
scale. This allows to isolate the presently uncalculable infrared dynamics in a set of 
universal parameters or functions. One is then led to the conclusion that relations of 
physical quantities can be calculated in perturbation theory. For certain observables, 
such as those derived from totally inclusive processes with no hadrons in the initial state, 
one obtains parameter-free predictions (except for the strong coupling) to logarithmic 
accuracy in the hard scale. 

Much effort has been invested in the computation of higher order QCD radiative 
corrections. For a few observables for which quark masses are unimportant, second 
order corrections in the strong coupling a s are known. Third order expressions exist for 
hadroproduction in e + e~-annihilation and r-lepton decays and for the Gross-Llewellyn- 
Smith and Bjorken sum rule in deep inelastic scattering. On the other hand, observables 
involving quark masses, for example decay widths of hadrons containing a heavy quark, 
are known only to first order. In most cases, the extension of present results is not 
simply a question of algebraic complexity and algorithms for handling it. The existing 
calculations have exploited presently known techniques to the frontier beyond which 
new methods need to be designed, a task which might not be completed soon. In this 
situation it is interesting to explore possible sources of systematically large corrections 
in higher orders, which, once identified, might then be taken into account exactly to all 
orders. In doing so, one may hope to reduce the uncertainty due to ignorance of exact 
higher order coefficients. 

Resummations of this type are familiar and necessary for many problems involving 
disparate mass scales. Large perturbative corrections in higher orders are associated 
with a logarithm of the ratio of these scales and can often be summed with renormaliza- 
tion group techniques, which allow to obtain the dominant power of that logarithm to 
all orders by one-loop calculations. Although in practice it is not always clear whether 
the logarithms dominate the constant pieces - in particular, since the coefficients of log- 
arithms grow geometrically with order of a s , whereas the constants grow factorially -, 
such a resummation is controlled by an "external" parameter (the logarithm of two 
scales) that can be varied, at least in fictitious limits. For the problem at hand, we as- 
sume that such renormalization group improvement has already been done or consider 
observables that depend only on a single scale. We are thus interested in systematically 
large contributions to constant terms with no external scale at our disposal. 

The estimation of constant terms of yet uncalculated coefficients, which is closely 
related to the problem of scale- and scheme-setting for truncated perturbative expan- 
sions, is often thought to attempt the impossible. While indeed for any particular 
observable it is impossible to assess with rigour the quality of a certain scale-setting 
prescription, some prescriptions may be supported by general physical arguments and 
turn out to be closer to exact results a posteriori. The most distinguished prescription 
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Figure 1: A typical diagram with multiple fermion loop insertions into the lowest order 
correction to a generic physical quantity. 

of this kind has been formulated by Brodsky, Lepage and Mackenzie 0. Observing 
that in QED all scale-dependence of the coupling results from photon vacuum polar- 
ization, they suggested that the effect of fermion loop insertions into a photon line in 
higher orders be absorbed into the scale of the coupling of a given order and proposed 
to apply this criterion to QCD as well. In practical applications, this suggestion has 
mainly been realized in second order in a s . To this accuracy, only a single fermion 
loop insertion needs to be calculated and the relevant contribution can be traced by its 
dependence on the number of light fermions Nf. In general, at order a™ +1 , the effect of 
one-loop evolution of the coupling can still be obtained by the highest power of Nf in 
the flavour-dependence of coefficients. The highest power of Nf originates solely from 
diagrams with n fermion loop insertions, such as in Fig. |1|, which are much easier to 
calculate than flavour-independent terms. In QCD, the identification of contributions 
related to one-loop renormalization of the coupling implies that the coefficient of the 
highest power in Nf should in fact not multiply Nf to some power, but the combination 
Nf — 33/2, as it appears in the Gell-Mann-Low function to leading order. Computing 
diagrams as in Fig. [I] and then replacing Nf by Nf — 33/2, one takes into account par- 
tial contributions from other diagrams, which are much harder to evaluate exactly (the 
precise identification of these contributions with particular diagrams is not straight- 
forward and gauge- dependent). It turns out that in all cases where comparison with 
exact second order results is possible, this replacement approximates the exact coeffi- 
cient amazingly well in the MS-scheme ||. In what follows we shall refer as "Naive 
Nonabelianization" (NNA) || to the hypothesis that a substantial part of higher order 
radiative corrections can be accounted for by running of the coupling, in the sense of 
substitution of Nf — > Nf — 33/2 in the term with the highest power of Nf. We should 
mention that BLM scale-setting and NNA are technically identical and the distinction 
is rather a matter of interpretation: The BLM scale-setting by its physical motivation 
does not necessarily claim that genuine higher order corrections missed by the above 
substitution should be small, whereas NNA assumes this stronger assertion. 

This paper is concerned with an exposition of details of the recent proposal M of 
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two of us to extend the BLM procedure to higher orders and, in particular, to sum 
all contributions associated with one-loop running of the coupling. We emphasize that 
this procedure has repeatedly been suggested in the past [|T], [| , but has apparently not 
been pursued to the point of practical implementation. A summation similar to, but 
technically different from Q was also presented in |J. We believe that this resum- 
mation is useful, because the effect of evolution of the coupling is a systematic source 
of potentially large perturbative coefficients. In such a situation it is advantageous 
to incorporate this effect in theoretical predictions even if it transcends a fixed-order 
perturbative approximation. 

To become definite, we use the following convention for the QCD /3-functionQ: 

P{a a ) = P a 2 s + (3ia 3 s + . . . ; 

-i( u -H> A =-(4^( io2 -f^)- (li> 

For a generic physical observable R (which for the sake of illustration we assume to 
depend only on a single, large scale Q), we may eliminate the iV/- dependence in radiative 
corrections in favour of the /^-function coefficients /3q, (3\ etc. For now (and most parts 
of this paper) we restrict ourselves to one-loop running, induced by (3$. The (truncated) 
perturbative expansion of R is written as 

N 

R-Rtree = I r n0 + r nX N f + . . . + r nn Nf]a s (Q) n+1 

n=0 

N 

= r a s (Q)J2[Sn + d n (-f3 ) n ]a s (Q) n , (1.2) 
n=0 

where 

d n = (-6vrr^. (1.3) 

Let us stress again that d n is unambiguously determined by diagrams with n insertions 
of fermion bubbles as in Fig. |l| and only S n requires a true n + 1-loop calculation. The 
5 n (n > 2) could be further separated into contributions from two-loop evolution and 
the effect of one-loop evolution on the genuine two-loop corrections (see for <5 2 ). We 
introduce 

N 

M N (-(3 a(Q)) = 1 + E d n (-(3 a(Q)) n (1.4) 

n=l 

^ince we discuss corrections proportional to (— f3oa s ) n , the reader accustomed to the normalization 
11 — (2/3)7Vj may think of an expansion in terms of a s /(4n). 
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as a measure of how much the lowest order correction is modified by summing N 
one-loop vacuum polarization insertions. To make an explicit connection to BLM scale- 
setting, we define 

a s (Q* N ) = <*,(Q) M N (-(3 a s (Q)) (1.5) 

absorbing the effect of vacuum polarization in the scale Q* N of the coupling. For N = oo, 
0£ a {Q*) (Q* = Qlo) is transparently interpreted as the running coupling averaged over 
the lowest order corrections. Isolating the integration over gluon virtuality in the lowest- 
order diagram, we may write 

r a s (Q) = a s {Q) J d 4 kF(k,Q) ^. (1.6) 



rn*\ r a s (Q) (a±ut?(i n N a s (fcexp[C/2]) 
r a s {Q) = - 2 ; ^ 9N = / d kF{k,Q) — , (1.7) 



Then 



l-f3 a s (Q)HQ* 2 /Qi) J K ^' k 2 

where C is a the scheme-dependent subtraction constant for the fermion loop. The 
difference between the leading-order BLM scale Q\ and Q* is precisely the difference 
between averaging the running coupling itself, or the logarithm of its argument, over 
the lowest-order diagram. The values of Q* and Q* N for iV > 1 depend on the value of 
a s (Q). Such a dependence is intrinsic to any generalization of the BLM prescription 
beyond leading-order and has previously been noted in R [7J. As a result we may 
express any physical quantity as 

R- Rtree = r (a(Q)M 00 (-foots (Q) ) + £ 5 n a s (Q) n . (1.8) 

Separating radiative corrections in this way and summing a partial set of them into 
Moo(— j3 a s ) can be motivated in different ways. First of all, as already mentioned, 
this procedure is supported empirically at the two-loop level. Provided the three-loop 
coefficient is sizeable, one might again expect that a significant portion of it can already 
be accounted for by M 2 (— /3 a s ) and, in higher orders, by M N (—j3 a s ). 

It must be noted that such a statement depends on the choice of scheme and, 
since the 5 n depend on Nf for n > 1, also on Nf for a given scheme. Although Q* 
is defined such that a s (Q*) is scheme- independent if one consistently uses the one- 
loop /9-function, an arbitrary ^/-independent finite renormalization can always ruin 
the supposed dominance of (3q a™ +1 -terms by implying a definition of the coupling that 
leads to anomalously large values of the remaining terms 5 n . This remark applies 
of course to any partial resummation, whether it is based on a systematic (external) 
parameter or not as in the case at handf]. In this respect a resummation utilizing one- 

2 We do not consider (3q on the same footing as, say, lnQ 2 /m 2 , since it is an intrinsic parameter 
of QCD. However, we shall somewhat loosely refer to the approximation of ignoring the S n as the 
"large-/?o limit". 
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loop running is similar to resummation of 7r 2 's that arise from analytic continuation 
from the space-like to the time-like region of, for instance, the Sudakov form factor |J 
or the coupling as in e + e~-annihilation or r decays ||KJ|| , and, in fact, comprises the 
resummation for r decays as a special case. Since empirically the /3o-term dominates 
second order radiative corrections for many observables in the MS scheme, one would 
not expect the resummation of /3oa™ +1 -terms to be numerically useful in schemes that 
differ from MS by a redefinition of the coupling that changes significantly the relative 
importance of flavour- dependent and flavour-independent terms. 

The dominance of (3q a" +1 -terms is lended further support from the behaviour of 
perturbative coefficients in large orders (renormalons) . It is precisely the effect of 
vacuum polarization, which leads to the expectation that in sufficiently large orders, 
the coefficients should indeed have the form 

r n ~ K {afo) n n\n b . (1.9) 



The series is thus divergent and, provided it is asymptotic at all, cannot approximate 
the true result to any desired accuracy, even if the exact coefficients could be computed 
to arbitrary order. The evaluation of diagrams as in Fig. [1] provides some insight into 
this ultimate limitation of perturbative QCD and into the approach to the asymptotic 
limit. In this paper, however, we are not so much interested in large orders and effects 
suppressed by powers of the hard scale Q. Summation of effects associated with one- 
loop running of the coupling should prove useful in intermediate orders (say, n = 
2 — 6), when perturbation theory is still reliable. It is amusing to speculate whether 
the empirical dominance of (3q at the two-loop level indicates that one is already close 
to the asymptotic behaviour. The latter is derived from a saddle point expansion and 
the large contribution arises from internal momenta either much smaller or much larger 
than the external scale Q. Proximity to the asymptotic behaviour - if true - would 
indicate that the distribution of internal momenta is already well-approximated by a 
Gaussian even though the main contribution to the Feynman integral is not from small 
or large momenta. 

The consideration of large orders provides also useful guidance to the limitations that 
arise from the restriction to one-loop evolution effects. As a matter of fact, diagrams 
associated with one-loop evolution do not provide the correct constants b and K in 
the large-n behaviour. At large n, the effect of two-loop evolution on a single gluon 
line and one-loop evolution on two gluon lines becomes equally important and should 
eventually be taken into account^. In practical applications, we will often find that 
the series has to be truncated due to its divergence at rather low n. Thus, many 
diagrams which are required to establish the formal large-n behaviour never become 
relevant. Moreover, provided 8\ is already not large, one should also expect that vacuum 
polarization insertions into two gluon lines will remain comparatively small. 

Let us repeat that the resummation of (3q a" +1 -terms can only be judged a posteriori 
and can fail to provide a good estimate of higher order corrections for any particular 

3 For example, an infinite number of two-loop insertions is necessary to obtain the correct value of 
b & 
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quantity. Even in this case, we believe that this resummation is physically motivated 
and that it is appropriate to absorb this particular class of higher order corrections into 
the scale of the coupling in lowest order. Application of BLM scale-setting to leading 
order (absorbing the contribution from one fermion loop) yields Q\ = Qe~ dl ^ 2 . Upon re- 
expansion of a s (Q*) to one-loop accuracy for the /3-function, this pretends a geometric 
growth d® LM = d™ to be compared with the factorial growth of the exact coefficients. 
Resummation to all orders corrects for this discrepancy by adjusting Q*^ such that the 
expansion of a s (Ql ) gives the correct values of d n . Thus to the accuracy of one-loop 
evolution the result of resummation is equivalently expressed as a s (Q)M oo (—/3 a s ) and 
we often prefer this form of presentation. It is no longer equivalent to a(Q^), if one 
adopts two-loop accuracy for a s . Since for large n the true d n will always outgrow d™ , one 
might conclude that the usual BLM scale-setting underestimates higher order radiative 
corrections. In practice, the effect is often just the opposite: The most important 
contributions come from intermediate n and it turns out that in many cases d\ is 
comparatively large {Q\ is small), so that d™ > d n in this region. The usual BLM scale- 
setting therefore typically overestimates the size of radiative corrections associated with 
one-loop running. 

The remainder of the paper is organized as follows: In Sect. 2 we develop in detail the 
technique to implement the resummation. We use a dispersion technique to reduce the 
problem to a one- dimensional integral over lowest order radiative corrections computed 
with finite gluon mass which is suitable to numerical evaluation. Thus compared to 
the complications of a complete higher order calculation, this resummation can be 
performed with little computational expense. Compared to a similar implementation 
of the standard BLM scale-setting fl2 |, the computation of higher order coefficients 



d n comes with no additional expense at all. It is convenient to introduce the Borel 
transform as a generating function of higher order radiative corrections. The principal 
value definition of the Borel integral serves as a starting point for the summation of the 
series within a certain accuracy, indicated by the presence of infrared renormalons. We 
shall also see that this definition requires all kinematic invariants to be large compared 
to Aqcd, reflecting the inapplicability of perturbation theory as a starting point for 
summation, if this requirement is not met. In this Section we also generalize the results 
of to quantities with anomalous dimension and include quarks with finite masses in 
the loops. 

In Sect. 3 we apply the resummation to the hadronic width of the r lepton and 
find a 10% decrease of a s {m T ) due to four- and higher loop corrections. We discuss the 
possibility of 1/m^-corrections in the light of our resummation. In general we prefer 
to be agnostic about power corrections and stick to perturbation theory. An exception 
to this rule is that we do not want to introduce power corrections in conflict with 
the operator product expansion in euclidian space without good reasons (which we do 
not have). Under this assumption we show that principal value resummation does not 
introduce 1/m^-terms to the decay width, provided the coupling is chosen appropriately. 
Sect. 4 contains a detailed discussion of the resummation of (— /3o) n a™ +1 -terms for the 
pole mass of a heavy quark. We combine the resummation with the exact two-loop result 
to give an estimate of the difference between the pole mass and the MS-renormalized 
mass. We keep finite quark masses inside loops which allows to trace the origin of large 
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coefficients to the relevant regions of internal momentum. 

In Sect. 5 we formulate one possible extension of our resummation to incorporate 
partially the effect of two-loop running resulting in (3\[3q ~ 2 -corrections. This extension 
is again guided by the flavour-dependence of coefficients and can be considered as an 
extension of recent work by Brodsky and Lu 0. The size of corrections is illustrated 
for the vector correlation function relevant to r decays and for the pole mass of a heavy 
quark. Sect. 6 contains conclusions. Three Appendices deal with technical issues: In 
Appendix A we derive a simple expression for subtractions needed for logarithmically 
ultraviolet divergent quantities and collect the details of results presented in Sect. 2. 
Appendix B contains analytic formulae for the lowest order radiative corrections to 
R e + e - and hadronic r decays with finite gluon mass. In Appendix C we list exact 
results for some abelian five-loop diagrams to the vector correlation function. 

In a companion paper Jl3| we shall discuss the implications of resummation for 
heavy quark decays and the determination of |V& C |- 

2 Techniques 

The aim of this Section is to develop a systematic approach to the calculation of dia- 
grams with an arbitrary number of fermion loop insertions, such as in Fig. 1. We assume 
a generic physical (short-distance) quantity R and a renormalization scheme that does 
not introduce an artificial Nf dependence, such as MS. It is also assumed that lowest 
order radiative corrections to R do not involve the gluon self-coupling. Subtracting the 
tree-level contribution, we are left with the (truncated) perturbative expansion 

N 

R-Rtree = Y,r n a^ 1 . (2.1) 

n=0 

The coefficients r n are polynomials in Nf , 

r n = r n0 + r nl N f + ...+ r nn N r f, (2.2) 

and we will calculate the coefficients r nn , which originate from n fermion loop insertions 
into the lowest-order diagram. We then write 

Tn = r [5 n + (-poTdn] , (2.3) 

where d n = (—Qn) n r nn /r absorbs the term with the largest power of Nf. The effect of 
one-loop evolution of the coupling on lowest order radiative corrections is then entirely 
contained in the d n and in the remainder of the paper we will consider the 8 n as correc- 
tions to the approximation of "Naive Nonabelianization" . As a measure of how much 
the lowest order radiative correction is modified by including N vacuum polarization 
insertions, we define 
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N 



M N (-(3 a s ) = l + ^d n (-/? a s ) n , 



n=l 



M oo (-/3 a s ) = Mjv^ooC-Aa,)- (2.4) 



Taken literally, the limit N — > oo does not exist, because the series diverges. It will be 
interpreted in the sense of Eq. ( |2.8|) below. For further use, we introduce the shorthand 
notation 

a s (fi) = -Po<x 8 (n) , (2.5) 

where \x is the normalization scale, which often will be suppressed for brevity. In this 
Section we are concerned only with technical aspects of the calculation of M^ia^) and 
Moo (a,). 

2.1 Borel transform vs. finite gluon mass 

A convenient way to deal with diagrams with multiple loop insertions is to calculate 
the Borel transform 



which can be used as a generating function for the fixed-order coefficients [fL4|: 



rn = (-Po) n £^B[R](u)^ . (2.7) 

Another advantage of presenting the results in form of the Borel transform is that the 
result for the sum of all diagrams can easily be recovered by the integral representation 



1 f°° i 

r M oo (a s ) = - due- u l a °B[R](u), 2i 

8s JO 



where the integration goes over positive values of the Borel parameter u. Note that we 
define the Borel parameter u with an additional factor (—/3q) compared to the conven- 
tional definition. In fact, Eq. ( |2.8|) requires some explanation: As it stands the integral 
is not defined, because the Borel transform generally has singularities on the integration 
path, known as infrared renormalons. We shall adopt a definition of the integral based 
on deforming the contour above or below the real axis or on a principal value prescrip- 
tion. These prescriptions are not unique and their difference, which is exponentially 
small in the coupling, must be considered as an uncertainty which can not be removed 
within perturbation theory. This will be discussed in more detail. A second question 
concerns the existence of the principal value integral and the behaviour of the Borel 
transform at u = oo. If we consider a physical quantity that depends only on a single 
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scale Q, then, to one-loop running accuracy, renormalization group invariance entails 
that the Borel transform can be written as {^ 2 /Q 2 ) u times a Q- and /i-independent 
function F{u). Combining the factor (fi 2 /Q 2 ) u with e~ u ^ as ^ in Eq. ( |2.8| ), we deduce 
that the principal value integral exists, provided that Q 2 is sufficiently large compared 
to Aq CD and F{u) does not increase faster than any exponential as u approaches posi- 
tive infinity Since the second property is satisfied in all examples which we shall meet 
(and we may conjecture that this is generally true for the Borel transform computed 
from higher order corrections due to vacuum polarization), we shall assume in general 
that all kinematic invariants on which R depends explicitly are sufficiently large com- 
pared to Aq CD (to be precise, the difference of two invariants is also an invariant). This 
is not an additional assumption needed to take the Borel integral. Without it pertur- 
bative methods are not applicable and there is no perturbation theory to start with. In 
particular, we do not know whether the Borel transform, which is useful in connection 
with short- distance expansions, can serve as a bona fide starting point to summation in 
the strong coupling regime^]. 

In simple cases the Borel transform can be calculated directly. This is due to the fact 
that the evaluation of diagrams with multiple fermion bubble insertions in Landau gauge 
corresponds to the evaluation of the lowest-order diagram with the effective propagator 

where 

H(k 2 ) =a s ln( -^e c ) (2.10) 



and C is a scheme- dependent finite renormalization constant. In the MS-scheme C = 
-5/3, in the V-scheme 0(7 = 0. 

For a chain of fermion loops contributing to a physical amplidude with euclidian 
external momenta, one can separate the integration over the gluon momenta to write 
it as 

roa^MM = J ^^Q) ^^ , (2-H) 

where the transverse projector that appears in the gluon propagator in Landau gauge 
is assumed to be included in the function F(k, Q) and Q stands for a collection of 
kinematic invariants. 

A crucial simplification arises since for diagrams with only one fermion bubble chain 
the Borel transformation applies to the expansion in a s of the propagator in Eq. ( |2.9| ) 
itself, rather than to the set of diagrams as a whole. The effective (Borel-transformed) 



gluon propagator is |14 



If some kinematic invariants are small, one might still be able to define the Borel integral as an 
analytic continuation. However, in this regime all power corrections in a short-distance expansion are 
of the same importance and the analytic continuation is useless unless the summation of the short- 
distance expansion is understood. 
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B[a s D™{k)]{u)=i5«"[p j " k f_~^r • ( 2 - 12 ) 



->AB/r„\i/„.\ _ - X AB 



Thus, the task of calculating the Borel transform of Feynman diagrams with bubble 
insertions reduces to the calculation of the leading-order diagram with the usual gluon 
propagator raised to an arbitrary power, which is familiar from analytic regularization 
r5|. This trick suffices to derive an expression for the Borel transform of the polarization 



operator with light quarks Jl4|, [16|, [17j , of the heavy quark self-energy |T7| and several 



more complicated cases in connection with heavy quark expansions, which can be found 

in ii7i, s a ■ 

However, in most phenomenologically interesting cases an analytic expression for the 
Borel transform is difficult to obtain, especially for observables involving more than one 
scale. Even if the exact Borel transform is obtainable, taking n derivatives to evaluate 
fixed-order coefficients (see Eq. ( |2.7| )) may turn out to be a complicated task. In this 
paper we shall work out a different technique, which extracts the desired information 
on higher orders from the lowest-order diagrams, calculated with a finite gluon mass 
jjj [12|]. Calculations of one- loop diagrams with finite gauge boson mass have become 
routine for electroweak radiative corrections, which allows to hope that this technique 
is generally applicable to a wide range of observables in QCD. In this way one obtains 
a concise integral representation for the Borel transform. 

For a while let us restrict our discussion to euclidian quantities, which are not 
sensitive to the gluon self-coupling to leading order .0 Call r (A 2 ) the leading-order 
radiative correction calculated with finite gluon mass A and r = r (0). To be precise, 
we define ro(A 2 ) as the sum of all Feynman diagrams to leading order, which in general 
may be ultraviolet divergent and need to be renormalized. In this Section we restrict our 
discussion to cases where no explicit renormalization is needed, which is the case, e.g. for 
the derivative of the polarization operator in Eq. (|3.1| ), or transition amplitudes related 
to heavy quark decays with on-shell mass renormalization. It is easy to show that this 
assumption is equivalent to the requirement that r (A 2 ) vanishes as A 2 — ► oo. Equally, in 
the case of Borel representation for the diagrams with fermion bubbles, we assume that 
fermion loops are renormalized, and no additional explicit renormalization is necessary. 
This assumption can easily be relaxed. A detailed discussion of renormalization is given 
in Appendix A, where we work out the missing overall subtractions for the general case. 

We keep the standard gauge-fixing and work with the propagator 



i5 AB ' 



k 2 - \ 2 + 



ic 



9iw ~ (1-0- 



k 2 -t\ 2 + 



(2.13) 



The relation to the Borel transform of the massless propagator in Eq. fl2.12| ) (in Lan- 
dau gauge, £ = 0, or Feynman gauge, £ = 1) is established by an (inverse) Mellin 
representation 

5 With the latter restriction, we do not have problems with gauge invariance, which otherwise 
prohibits introduction of a finite gluon mass, unless QCD is embedded, for example, in an SU(3) 
Higgs model. 
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1 



1 1 



-1/2+ioo 



k 2 - X 2 2m k 2 



J duT(-u)T(l + u) I- 



A^" 
k 2 



(2.14) 



-1/2-ioo 



Writing down the leading-order radiative correction with non-vanishing gluon mass A 
(in Landau gauge) as 



(2.15) 



ro(A^) = / d*kF(k,Q)- 



k 2 - X 2 

and comparing to Eq. Q2.11|) , one finds the identity |2C 

-1/2+ioo 



ro(A) 



— J duT{-u)T{l+u)[-^e c \ B[R](u) (2.16) 



-1/2-ioo 



Taking the inverse, we get [§J 



B[R](u) 



sm(iru) 

7T 

sm(iru) 
nu Jo 



00 dX 2 



'A 2 



r (A 2 )-r (0) 



dX 2 r' (X 



(2.17) 



where r' (X 2 ) = (d/dX 2 )r (X 2 ). If r (A 2 ) ~ A 2a for small A 2 , the first line exists, when 
< u < a. The second line provides the analytic continuation to an interval about 
u = 0, so that derivatives at zero can be taken. Thus fixed-order coefficients r n can be 
expressed in terms of the integrals 



J k = dX 2 ln k (X 2 /fi 2 )r' (X 2 ) k<n. 
Jo 



For instance (cf. Eqs. (|2l| and (p??))) 



do 
di 

d 2 
d 3 



- l-M = 1 

ro 

- [Ji + C J ] 

^0 

1 

ro 
1 

ro 



-Jo - 2C,h 



7T 



J 3 + 3CJ 2 + (3C 2 - vr 2 ) J x + (C 3 - Ctt 2 ) J 



(2.18) 



(2.19) 



etc. and r n = ro (—(3o) n d n . r\ coincides with the result of Smith and Voloshin |L2] (who 
use a scheme with C — 0), integrating of their expression by parts. 
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Next, we evaluate the Borel integral in Eq. ( |2.8|) to obtain a compact answer for the 
sum of diagrams with fermion bubble insertions to arbitrary order, up to the ambiguities 
caused by renormalons as mentioned below Eq. fl2.8Q . For the following derivation we 
assume that the Borel integral is defined with a contour slightly above the positive 
real axis. We also recall that this integral exists, if all kinematic invariants are large 
compared to Aq CD , which we assume. To find a representation for the so-defined sum 
in terms of an integral over A 2 , one needs to insert Eq. ( |2.17| ) in Eq. ( |2.8| ) and take 
the li-integral explicitly. The u-integration is elementary, but the interchange of orders 
of integration in u and A 2 can not be done immediately, because the A 2 -integral in 



Eq. ( p,17|) is not defined for all u on the integration path. It is convenient to use the 
first representation in Eq. ( p. 17 ) and write it asf] 



r M oo (a s ) 



+ 



1 

2ni Jo 



1 

2tH 



due~ u/as 



u/a s 


f 00 d\ 2 ( 




/o "A 2 " 1, 


dX 2 




A 2 





/i- 

— u 



r-o(A 2 ) - r (0) 



r (A 2 )-r (0) 



(2.20) 



The integration over u in the first of the two terms above can easily be taken, rotating 
the contour to the positive imaginary axis. The quarter-circle at infinity does not 
contribute, because of the behaviour of the Borel-transform at infinity as discussed 
above. The A 2 -integral exists for positive imaginary u and the order of integrations 
can be interchanged. The u-integral in the second term would equally easily be taken 
by rotating the contour to the negative imaginary axis, but the infrared renormalon 
singularities on the real axis obstruct this deformation. It can only be done at the price 
of picking up the residues from the infinite number of poles. A more elegant solution 
is to rotate first the integration contour over A 2 to the second sheet: A 2 — > X 2 e~ t(jT+t \ 
Here we note that for euclidian quantities, there are no singularities of the A 2 -integrand 
in the lower complex plane. Next, the ^-integration can again be performed by rotation 
to the positive imaginary axis, as above. This integration gives 



MOO 

/ du exp{-u[l/a s + ln(|A 2 |//i 2 e c )] - ie} = — 
Jo 1 



+ a s ln(|A 2 |/7i 2 e c ) - ie 
and introduces a pole singularity, located at 

A| = -fJ 2 exp[-l/a s - C] , 



(2.21) 



(2.22) 



which is simply the position of the Landau pole in the running coupling (in the V- 
scheme). Now one can rotate the A 2 -integral back from the second sheet to the real 

6 The separation of the two terms in the following equation reintroduces the pole at u = in each of 
the two A 2 -integrals. The pole cancels in the sum of both terms, and both terms can be manipulated 
separately regardless of this pole. This can be justified as follows: One splits the A 2 -integral before 
separation into the two terms at some sufficiently large X^. The contribution from A^ to infinity can 
be handled without difficulty. For the integral form to A 2 - one can proceed as following Eq. ( 2.2(.l| ) 
and both terms have no singularity at u = 0. 
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Figure 2: One-loop running coupling in the V-scheme, C — 0, (broken line) and 
effective coupling — $(A 2 ) (solid line) as functions of A 2 /Ay . 



positive axis. In this way one encounters the pole in Eq. fl2.21|) , whose residue has to 
be added. Collecting all the terms, we get after some algebra: 



r M oo (a s ) 



dX 2 



A 2 |l + a s ln(-A 2 /^ 2 e c )| 2 



r (A 2 )-r (0) 



1 

+ — 



r o (A|-k)-r o (0) 



(2.23) 



where the ie prescription corresponds to defining the Borel integral above the real axis: 
The contours are rotated in the opposite directions if the Borel integral is defined below 
the real axis, with the only modification that the sign of the ze-prescription is reversed 
in the result. Finally, integrating by parts in the first term, we obtain 

WaOMooM//)) = / dA 2 $(A 2 K(A 2 ) + [r (A!-k)-r (0)] (2.24) 
jo 



where 



<£(A ) = arctan 

7T 



a s (yLt)7T 



1 + a.Ox) ln(A 2 /// 2 e 



2 e C\ 



e(-A| - a 2 ) 



(2.25) 



which coincides with the result of 0], obtained by a different method. Note that the 
term with the B-function exactly cancels the jump of the arctan at A 2 = — A 2 . 

Eq. ( p.24|) presents the desired answer for the sum of diagrams with any number of 
fermion bubbles in terms of an integral over the gluon mass. This relation is one of 
the main technical tools which we suggest in this paper, and we discuss its structure in 
detail. 
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First, notice that Eq. ( |2.24|) has a very intuitive and simple interpretation: The 
quantity ?"q(A 2 ) (with certain reservations) can be considered as the contribution to the 
integral from gluons of virtuality of order A 2 , and the function — $(A 2 ) can be understood 
as an effective charge. At large scales, — <3>(A 2 ) essentially coincides with ov(A), the QCD 
coupling in the V-scheme [[[]], but in distinction to it remains finite at small A 2 , see Fig. 
The absence of a Landau pole in this effective coupling implies that the integral is a 
well-defined number and the fact that we have started with the expression in Eq. ( |2.8| ) 
that is ill-defined due to infrared renormalon singularities ( equivalent ly, attempted to 
sum a non-Borel summable series) is isolated in the Landau pole contribution r (A|). 
Whenever infrared renormalons are present, r (A 2 ) develops a cut at negative A 2 . The 
real part of the above expression for the bubble sum coincides with the principal value 
of the Borel integral and, in particular, coincides with the Borel sum, when it exists. 
The imaginary part provided by r (\ 2 L ) coincides with the the imaginary part of the 
Borel integral, when the contour is deformed above (or below) the positive real axis. We 
therefore adopt the real part of Eq. (|2.24j) as a definition of Moo(a s ), and consider the 
imaginary part (divided by 7r) as an estimate of the intrinsic uncertainty from summing 
a (non-Borel summable) divergent series without any additional nonperturbative input. 
Note that this imaginary part is proportional to a power of X 2 L /Q 2 ~ exp(— l/a a (Q)) 
and therefore is suppressed by powers of Aqcd/Q- 

Second, notice that the result in Eq. fl2.24j ) is manifestly scale- and scheme-invariant, 
provided the running of the coupling is consistently implemented with the one-loop (3- 
function: a s (/xi) = a s (//2)/(l + ^(^2) ^(^1/^2))- I* 1 particular, a s {Q*) is formally 
independent of the finite renormalization C for the fermion loop. To be precise, in 
schemes that do not introduce artificial flavour dependence, the coefficients of the ex- 
pansion that relates the couplings in two different schemes have itself an expansion of 
the type of Eq. (|2.2|) . Within the restriction to fermion loops with subsequent restora- 
tion of /3q one must again keep only the highest power of Nf in these coefficients. If 
couplings are related in this way, the numerical value for a s {Q*) = a s (Q)M 00 (a s (Q)) is 
scheme- independent. 

In practice, working with the full QCD /^-function one introduces a scheme- 
dependence in higher orders in 1/Nf similarly as usually with higher orders in a s (Q). 
The quality of the NNA as compared to exact coefficients depends on the numerical 
magnitude of terms that are formally suppressed by powers of Nf (or (3 ) and there- 
fore depends on the scheme. Since empirically the success of NNA in low orders is 
observed in the MS scheme, we cannot expect it to hold in schemes that differ from 
MS by terms that are formally of higher order but numerically large. Thus, although 
the resummation accounts exactly for the contribution of fermion loops in any scheme, 
we believe that its phenomenological relevance is tied to the MS scheme, or others that 
are "reasonable" in the above sense. 

Third, we remark that Eq. ( |2.24|) applies without modification to quantities like in- 
clusive decay rates, which can be obtained starting from a suitable amplitude in euclid- 
ian space and taking the total imaginary part upon analytic continuation to Minkowski 
space. The structure of the A 2 -integral remains unaffected, and it is only the quantity 
r(,(A 2 ) that should be substituted by the corresponding decay rate calculated with finite 
gluon mass (For heavy quark decays, no explicit renormalization is needed, when the 
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decay rates are expressed in terms of pole masses.). Similarly, Eq. ( p. 24 ) is applicable to 
various semi-inclusive quantities, with the only formal restriction that a weight-function 
chosen to specify the final state does not resolve quark- ant iquark pairs in fermion bub- 
bles, that is, their phase space must be completely integrated. It is not applicable to 
quantities like hadronic event shapes. 

The derivation of Eq. ( 2.24 ) is only slightly modified, when ?"o(A 2 ) represents a 
physical cross section. In this case, r (A 2 ) is written as a sum of virtual and real gluon 
emission, where real gluon emission occurs only when A 2 is below a certain threshold 
At(Q) 2 , 



r (A 2 ) = rw(A 2 ) + r real {\ 2 ) 6(A T (Q) 2 - A 2 ) 



(2.26) 



When all kinematic invariants, represented by Q, are large compared to Aq CD , the same 
is true for Xt{Q) 2 - It is simple to show that Eq. ( 2.1 7|) holds unmodified for the Borel 



transform of a cross section with Eq. (|2.26| ). When deriving Eq. (|2.24| ) one splits the 
A 2 -integral in Eq. ( [2.17| ) at \t(Q) 2 - The contribution from to \ T {Q) 2 can be treated 



in the same way as before (adding and subtracting a contribution from a semicircle with 
radius Xt{Q) 2 )- For the contribution from Xt{Q) 2 to infinity, one may interchange the 
order of A 2 - and u-integrals directly, since the first exists for all u. Then, the ^-integral 
can be taken, since Xt(Q) 2 is large compared to Aq CD (and therefore larger than — \ 2 L ). 
After the w-integral is taken, we combine both pieces. The final result is identical to 
Eq. flpD. 

We should also mention that the Borel transform of physical cross sections gener- 
ically contains a factor sin7ra, which seems to obstruct the rotation of the w-integral 
to the imaginary axis. Closer inspection of Eq. ( |2.20| ) shows that the two separate 



A -integrals there correspond in fact to the Borel transform of 



B[R]{U) , (2.27) 
sin iru 



which is sufficient to guarantee convergence at the imaginary u relevant to the two 
terms in Eq. ( |2.20| ). 



2.2 Quantities involving renormalization 

So far the derivation was restricted to quantities, where the lowest order radiative 
correction r$ is ultraviolet (and infrared) finite. Then, in diagrams involving insertion 
of fermion loops into the gluon line, counterterms are needed only for the fermion loop 
subdiagrams, but no overall subtraction for the diagram itself (after summation of all 
diagrams that contribute to r n ) is required. We consider now the more general case, 
that the quantity R(Q) has an anomalous dimension. One can still use Eq. ( |2.9| ), but 
the resulting Borel transform is singular at u — 0. This singularity is compensated by 
adding a function Sr(u), which accounts for the missing counterterms. The function 
Sr(u) in the MS scheme is given in Eq. (|A.7| ). 
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When expressed in terms of the lowest order radiative correction with non-zero 
gluon mass, the necessity of additional subtractions is reflected in a divergence of the 



A 2 -integral at large A 2 for u = in the second line of Eq. ( |2.17 ), since r (A 2 ) grows 



logarithmically at large A 2 (We assume a logarithmic ultraviolet divergence.). The 
additional counterterms not associated with fermion loop insertions amount to the 
subtraction of the leading term in the large-A 2 expansion of r (A 2 ) and to the addition 
of some scheme specific contributions. In the MS scheme, Eqs. ( |2.17| ) and ( |2.24j ) are 
replaced by (C = —5/3) 



B[R](u) 



smnu i ,.o 
dX 

TTU 



- G (u) 
u 



sin TTU 



TTU 



A 2 



e(A 2 



(2.28) 



and 



roasM^as) = ^ d A 2 $(A 2 ) ^(A 2 )-^6(A 



2 ^ e - c : 



+ / — (Go(«)-roo)+r c 
J u 



arctan(7ra s ) 1 , / 9 9 

+ - In 1 + vr 2 a 2 

7ra« 2 V 



[ro(Al)-r (0)] 



(2.29) 



Here, as always, a s = a s (fi). The derivation of this result together with the definition 
of all new quantities that appear in the above equation is given in Appendix A. In 
particular, Gq(u) is essentially the anomalous dimension of R (to leading order in 1/Nf). 



2.3 Quarks with finite masses 

In practical applications it can be important to trace the number of active fermion 
flavours, which may effectively decrease in high orders because important integration 
regions shift towards decreasing momenta, if the series is dominated by infrared renor- 
malons in large or intermediate orders. Thus, it is worthwhile to generalize the above 
technique to include massive quarks in fermion loops. For definiteness, we shall consider 
here the case of one massive, and Nf — 1 exactly massless flavours. The generalization 
to several massive flavours is then obvious. (3q is to be taken with Nf flavours, including 
the massive one, since we have in mind a situation, where the quark mass is non-zero, 
but smaller than the external momenta Q. We study the decoupling of quarks with 
finite masses in higher orders on a particular example in Sect. 4. 

The vacuum polarization IT(A; 2 ) in Eq. ( |2.9| ) includes summation over flavours. With 
one massive quark of mass m, it is modified to 



n(A; 2 ) = a s ln(-fe 2 //i 2 ) + C - A(k 2 ,m 2 ) 



Mk - m ) = <w%/ jrpW'.'"')- 1 ]. < 2 ' 3 °> 



16 



where 



p(s,m 2 ) : 

Performing the integral gives 



1 + 



2m 2 



\-—Q(s-Am 2 ) 



(2.31) 



A(k 2 ,m 2 ' 



1 



In 



+ 



In 



+ In 



Am 2 

1 

2 



4m 2 / 2m 2 
1 + 



A- 2 



A: 2 



2m' 



+ \ 1 



4m 2 



(2.32) 



To derive the expression for the sum of fermion loop insertions, we follow the method 
of ||, [L2| and substitute the effective propagator (1 + n(A; 2 ))~ 1 in Eq. fl2.11|) by the 
dispersion relation 



dX 



2 1 Imn(A 2 ) 



+ 



l + n(A: 2 ) Wo k 2 - A 2 |l + n(A 2 )| 2 k 2 - X 2 L W(X 2 L ) 
where IT(A 2 ) = d/d\ 2 Il(X 2 ) and X 2 L < is the solution of 

l + n(A|) = 



(2.33) 



(2.34) 



provided it exists. If no solution exists, the second term in Eq. ( |2.33| ) is absent. We 
now use 



1 



7T 



Imn(A 2 



l + ^-(p(A 2 ,m 2 )-l) 



(2.35) 



and write 



1 



1 / k 



1 • 



A; 2 - A 2 A 2 ^A; 2 -A 2 
Interchanging the order of integrations in k and A 2 , we arrive at 



(2.36) 



r Moo(a s ) = -a s 



dX 2 r (A 2 )-r (0) 

a 2 |i + n(A 2 )| 2 



1 



1 + — WV,mVl) 



r (A 2 )-r (0) 

+ A|n' ( A 2 ) (2 - 37) 



which coincides with Eq. ( |2.23 ) in the limit m — ► 0. 



A representation of the Borel transform which allows the calculation of fixed-order 
perturbative coefficients is obtained in a similar way. Starting from Eq. ( |2.11| ), one only 
needs 
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ImB 



i + n(A 2 ) 



(u) = l me -^(\ 2 )/as = e -uReIl(\l)/a s gin \_ u l mU ^ 2 )/ 



(2.38) 



Then we proceed as in the massless case and get 



B[R](u) 



1 

vr Jo 
( A 



x 



dX 



2 r i r 

- r (A 2 ) - r (0) sinju" 



1 + 



67T/?o 



(p(A 2 ,m 2 )-l)]} 



— e c \ exp{ttReA(A 2 ,m 2 )} . 



(2.39) 



For m 2 = we recover the old result in Eq. fl2.17| ). To avoid possible bad behaviour at 
A 2 — > oo and to separate the mass dependence, we add and subtract the expression for 
m 2 = 0. Defining 



T(u, A 2 , m 2 ) = sin {wtt[i + ^-(p(A 2 ,m 2 ) - 1)] } f-^e J exp {uReA(A 2 , m 2 )} , (2.40) 



we obtain, finally 



B[R](u) 



roc ( \2 ^ 

dA 2 ^e c 



vr Jo 



dA^ 
A 2 



-o(A 2 ) 



r (A 2 ) - r (0) T(«, A 2 , m 2 ) - T(n, A 2 , 0) 



(2.41) 



The second line gives the correction due to finite quark masses inside loops. We derived 
this result for quantities that do not need renormalization. However, in the MS scheme 
subtractions are mass-independent. Therefore, only the first term in Eq. ( j2.41j ) is 
affected by subtractions, and in precisely the same way as for m 2 = (cf. Sect. 2.2). 

We should point out that Eq. ( |2.37| ) for the sum of fermion loops has been derived 
through a dispersion relation. For massless quarks, comparison of the derivation of 
Eq. ( gjp with the one in |§] shows that the result obtained by a dispersion relation 
coincides with the principal value prescription of the Borel integral. We do not offer 
such an equivalence for massive quarks and there is reason to doubt that it is true in 
minimal subtraction schemes for the fermion loops. To illustrate this, suppose there 
were only a single massive particle inside loops, which produces a negative contribution 
to the /3-function. For any finite value of the mass of this particle, the vacuum po- 
larization behaves as k 2 /m 2 at very small virtuality k. Therefore the factorial growth 
of coefficients from the infrared region of integration is eliminated at very large order. 
There are no infrared renormalon poles and therefore no ambiguity in the Borel trans- 
form (though it may be sharply peaked at those values, where singularities appear as 
m approaches zero). On the other hand 1/(1 + U(k 2 )) does have a pole in the euclidian, 
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when m is smaller than a critical value and U(k 2 ) is defined by minimal subtraction. 
This leads to an ambiguity in the second term in Eq. (|2.37|) . The critical value is of 
order Aqcd and the discrepancy will be noticeable only for such small quark masses. 
However, for quark masses of order Aqcd, the quark mass must be considered as an 
infrared regulator and one can no longer straightforwardly identify an infrared param- 
eter, such as the gluon condensate, with non-analytic terms in a finite gluon mass. As 
is well-known, the gluon condensate must be redefined to absorb the non-analyticities 
in light quark masses as well. This is indicated by the highly singular behaviour of the 
Borel transform as the quark mass goes to zero. In our practical application the men- 
tioned discrepancy is numerically irrelevant, and we do not pursue this point further. 
Notice that such a difficulty does not appear in finite order perturbative coefficients, 
derived from Eq. ( |2.41| ), and affects only the Landau pole term, in which we are mainly 
interested as a measure of intrinsic uncertainty. 



2.4 Renormalon ambiguities and extended Bloch-Nordsieck 
cancellations 

We want to emphasize the close relation of renormalon singularities to non-analytic 
terms in the expansion of lowest-order radiative corrections in powers of the gluon mass 



20| . A formal relation between singularities in the Borel plane and non-analytic terms 
in A 2 is established by Eq. ( |2.16|) . Each non-analytic term proportional to [\fX 2 ) 2n+1 



in the expansion of ro(A 2 ) at small A 2 is in one-to-one correspondence to a single- 
pole singularity of B[R](u) at positive half-integers u — n + 1/2. Each non-analytic 
term proportional to A 2n In A 2 corresponds to a single pole at positive integer u = n. 

Likewise, non-analytic terms in the expansion at large A 2 of type a/A^ ^ or A~ 2n In A 2 
correspond to a single-pole singularities of the Borel transform at negative u = — n — 1/2 
and u = —n, respectively. For double or higher poles, higher powers of In A 2 appear. 

We recall that through the presence of singularities for real positive values of the 
Borel parameter the perturbation series signals its deficiency: Explicit non-perturbative 
corrections must be added to make the full answer unambiguous. In fact, the infrared 
renormalon problem is just one manifestation of a generally accepted wisdom: Pertur- 
bative calculations are only reliable if the essential integration regions include momenta 
much larger than Aqcd- The above relation between infrared renormalons and non- 
analytic terms in the expansion in powers of the infrared regulator like the gluon mass 
is thus natural and expected. A comment is necessary, however, to explain why only 
non- analytic terms in the expansion at small A 2 are related to the infrared behaviour, 
and simple power-like terms, A 2ri , are not. 

A small gluon mass A 2 ~ Aq CD not only eliminates contributions of small momenta 
k 2 ~ Aq CD , but also modifies the gluon propagator at virtualities of order k 2 ~ Q 2 . In 
this region of momenta 1/(A: 2 — A 2 ) can be expanded and produces (infrared insensitive) 
corrections of the form X 2 /k 2 ~ X 2 /Q 2 . In perturbative calculations these terms are 
unimportant, since the corrections are suppressed by powers of Q 2 . Hence the common 
practice to use the finite gluon mass as an IR regulator in calculations of various QCD 
observables. The famous Bloch-Nordsieck cancellations guarantee that terms propor- 



19 



tional to In A 2 , which appear at intermediate steps of the calculation, cancel in final 
answers for inclusive quantities. Calculations aiming at power-like accuracy must trace 
accurately the fate of power-like terms in the gluon mass. Since analytic terms A 2ri 
come entirely from the expansion of the gluon propagator at large virtualities of order 
Q 2 , they disappear when the regulator is removed. Only non- analytic terms are related 
to the infrared behaviour, and indicate the failure of perturbation theory to account for 
this region properly. Thus, tracing non-analytic terms in the expansion of leading-order 
radiative correction at small gluon masses, one can trace the sensitivity of particular 
quantities to the IR region, and, in particular, judge upon existence of non-perturbative 
corrections suppressed by particular powers of Q 2 . 

Stated otherwise, the absence of particular power-suppressed corrections in physi- 
cal quantities can be understood as an extension of the Bloch-Nordsieck cancellations. 
For example, the fact that in Wilson's operator product expansion for the correlation 
function H(Q 2 ) the l/<5 2 -corrections are absent, corresponds in this language to the 
cancellation of corrections of order A 2 In A 2 in the correlation function. Note once more 
that analytic terms proportional to A 2 do not necessarily cancel, and in fact they are 
present in a quantity, closely related to U(Q 2 ), the r-lepton hadronic width, see be- 
low. In turn, the existence of a gluon condensate contribution, {G 2 )/Q 4 , implies that 
contributions of order A 4 In A 2 do not cancel. 



3 Hadronic r decays 

In this Section we apply the summation of one-loop running coupling effects developed 
above to observables related by analyticity to the correlation function of vector (axial- 
vector) currents 

n F (?) = (g^-g 2 v) n W 2 )=» / d 4 xe^(0|T{it( x )^( )}|0), Q 2 = -q 2 . (3.1) 

Quantities of prime physical interest related to H(Q 2 ) are the cross section of e+e~ 
annihilation 

W.) ^ r ''T-~ "t° n f = iari-nv(.) (3.2) 

(neglecting Z-boson exchange) and the r-lepton total hadronic width 



I — > u T + hadrons) 



r(r- 



(3.3) 



= 12vr 



ds 



m% 



mi 



1 + 2-£- ) Im n« (a) + Im n* ' (s) 



where s = q 2 and we introduced the decomposition 
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LL V/A 



(<?) = ( 



^ + <^n^(<z 2 



(3.4) 



for the vector and axial- vector correlation function and defined E[W(s) = Hy (s) + 
n?(s). For the purpose of our discussion we neglect the strange quark mass and omit 
the overall CKM factor |V ud | 2 + \V US \ 2 » 1. 

The exact (nonperturbative) correlation functions should be analytic in the complex 
s-plane cut along the positive axis. Exploiting this property, we may transform Eq. 
into 



R T = Qiri 



ds 



mt 



1 + 2- 



mt 



nW( s ) + n(°)( 4 



(3.5) 



|s|=m| 



The same analytic property holds to any finite order in perturbation theory in a s (/i) 
(although the discontinuity is arbitrarily wrong at small s). Eqs. Q3.3| ) and fl3.5| ) are 
equivalent if the correlation functions are substituted either by the exact values or 
finite order perturbative expansions. In addition, in perturbation theory, the vector 
and axial-vector contributions coincide. The equivalence of Eqs. (|3.3|) and (|3.5Q does 
not hold in renormalization group improved perturbation theory or after fermion loop 
summation. This will be discussed extensively below. 

In what follows we concentrate on r decays, which despite the low energy scale 
involved, are considered to provide one of the most reliable determinations of the QCD 
coupling 1 22]. The state-of-the-art perturbative calculations yield |23| in the MS scheme: 




+ 



'a s {m 2 T Y 



7T 



(6.3399 - 0.3792JV, 



(48.5832 - 7.8795iV7 + 0.1579iV 



f> 



+ 0(at). (3.6) 



The leading non-perturbative corrections due to contributions of local operators of di- 
mension 6 have been estimated [211 p2| and turn out to be small, below 1% compared 



to the tree-level unity in brackets above. Thus, the principal uncertainty of the deter- 
mination of a s from the r hadronic width comes from unknown higher-order terms in 
the perturbative expansion. 

The purpose of this Section is twofold: We analyze summation of (— /3oa s ) n per- 
turbative corrections and conclude - with certain caveats - that the cumulative effect 
of higher order corrections beyond order is somewhat larger than the exact order 
a^-correction. Resummation of vacuum polarization reduces the value of a s (m T ) ex- 
tracted from t decays by about 10%. Second, we endeavour to clarify and speculate 
on a conceptual point: Several authors 



M| [25|, 22 1 considered (with different conclu- 
sions) the possibility that power corrections proportional to 1/m 2 . may creep into R T 
from various sources (summation of large orders in perturbation theory - infrared and 
ultraviolet renormalons -, freezing of a physical coupling in the infrared, violations of 
duality), whose absence in the approach reviewed above is crucial to ascertain its power 
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for the determination of a s . We elaborate on two aspects, which are sometimes omitted 
from the discussion: The necessity to define a coupling parameter to power-like (in 
l/m T ) accuracy and the distinction between hypothetical 1/m^-corrections to the OPE 
of correlation functions in the euclidian and to the r decay width itself. 



3.1 Higher orders in (—flo) a s 

To start with, let us demonstrate that Naive Nonabelianization provides indeed an ex- 
cellent approximation in low orders. The NNA approximation is obtained from Eq. ( |3.6| ) 
keeping the term with the leading power of Nf only, and restoring the full (3q by the 
substitution Nf — > Nf — 33/2. The result is 



R 



NNA 



+ (6.2568 - 0.37927V,) 



7T 



7T 



+ 



a "( mT A (42.9883 - 5.2107AT/ + 0.1579JVJ} 



+ 0(aj). (3.7) 



The accuracy is impressive: For the practical case Nf = 3 NNA predicts a coefficient 
28.7773 for the (a s /7r) 3 correction, to be compared to the exact coefficient 26.3658. 
We should mention, however, that this coincidence is also slightly misleading, since a 
substantial part of the exact coefficient is from a combination of /3-function coefficients 
and lower order coefficients of the correlation function Tl(s) and results from contour 
integration in Eq. fl3.5|) . 

To estimate higher orders, we recall that for the polarization operator, Eq. ( |3.1| ), 
an analytic expression is available for the Borel transform of the sum of diagrams with 



fermion loop insertions [14], [16], |26], |17|] . To get rid of an (irrelevant) overall subtraction, 
it is convenient to take one derivative with respect to Q 2 and to consider 



D(Q 2 ) = 



■ dn(Q 2 
dQ 2 



(3.8) 



Then (the simple representation quoted here is adapted from [IB], see [[26J, [T 

' 12 \"« - ( _ l) k k 



B[D(Q 2 )](u) = ^ l^e c 



and 



E 



l-{l-uf^ 2 {k 2 -{l-uff 



(3.9) 



B[R T ](u) = 127r 2 B[D(Q 2 )](u)sm(7ru) 



1 2 2 1 

TTU 7r(l — u) 7r(3 — u) 7r(4 — u) 



(3.10) 



Taking derivatives of B[R T \(u) (see Eq. ( p.7|) ) it is easy to evaluate fixed-order perturba- 
tive coefficients. In Table 0, second to fourth column, we give values for the coefficients 
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n 


Jt.MS 
a n 


d T ' v 


m;> ms 


a n 


M r,LDP 
n 





1 


1 


1 


1 


1.329 


1 


2.2751 


0.6084 


1.521 


0.6918 


1.578 


2 


5.6848 


0.8788 


1.819 


3.1035 


1.855 


3 


13.754 


-0.3395 


1.984 


2.1800 


1.898 


4 


35.147 


3.7796 


2.081 


30.740 


2.027 


5 


84.407 


-14.680 


2.134 


-34.534 


2.000 


6 


248.83 


99.483 


2.170 


759.74 


2.094 


7 


525.38 


-664.00 


2.187 


-3691.4 


2.041 


8 


3036.0 


5400.06 


2.210 


42251 


2.042 



Table 1: Coefficients for n fermion loop insertions into one- loop radiative corrections for 
the r-lepton hadronic width and partial sums of the perturbation theory for a s (m T ) = 
0.32 [a s (m T ) = 0.229]. See text. 

d n defined in Eq. fl2j for n < 8 in the MS- and V-schemeQ (C = -5/3 and C = 0), 
and partial sums of the perturbative series M n (a s (m r )) in the MS-scheme, defined in 
Eq. (|2.4| ), for a^ s (m T ) = 0.32 |E| and taking Nf = 3 massless flavours. 

Truncating the expansion of R T at its minimal term, one deduces from the fourth 
column of Table |1] that the effect of summation of (— /3oa s (m T )) n -corrections increases 
the leading order correction by factor My' M (0.229) ~ 2.19, which is to be compared 

with 1.803, obtained from the exact coefficients up to order a^. The value of Mj' M (a s ) 
from the truncated series is in good agreement with the result of resummationP], 

M^ MS (0.229) = 2.233. (3.11) 



Taken at face value, higher order vacuum polarization effects lead to a significant in- 
crease of radiative corrections beyond the a ^-approximation. The cumulative effect 
amounts to somewhat more than the a^-correction itself. 

Before we discuss its impact on the determination of a s , we note that it has become 
customary not to use a fixed order approximation for R r , but to employ the approach of 



7 Note that both series are far from the expected asymptotic behaviour. There are two reasons 
for this: First, for the current correlation function, the formal large-n behaviour is dominated by 
ultraviolet renormalons, whose suppression in the MS scheme with C = —5/3 is significant in low 
orders. Therefore, one does not expect to see sign-alternating behaviour in low orders. Second, the 
contour integration rearranges coefficients and postpones the onset of the asymptotic regime, see below. 

8 We have obtained this value in three different ways: (1) Starting from Eq. (2.24) with the leading 
order corrections to hadronic r decays with finite gluon mass collected in Appendix B; (2) Taking the 
principal value Borel integral of Eq. (3.1C) and (3) Computing the principal value integral of Eq. ( |3.9| ) 
along a circle of radius in the complex s-plane and then taking the contour integral along the circle 
according to Eq. (3.5). The intrinsic uncertainty due to IR renormalons is small, ±0.02, since the 
leading pole at u = 2 disappears in R T in the large-/3o limit. 
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Le Diberder and Pich (LDP) [110], who resum exactly the effect of the running coupling 
along the circle in the complex plane, but not on Q 2 dU/dQ 2 itself. This resummation is 
motivated by the observation that evolution along the circle generates a series of large 
higher order corrections, which is convergent, but only barely so at the actual value of 
a s (m T ). The resummation of Le Diberder and Pich (restricted to one-loop running of 
the coupling) is automatically included in our resumation, since the running coupling 
is not expanded in a s (m T ) in the derivation of Eq. ( |2.24j ). It is yet instructive to apply 
the procedure of (TpJ to fixed order approximations for Q 2 dU/dQ 2 . The successive 
approximations for R T are then given by 



N 

M T / DP \a s {m T )) = Y, d n A n{a s {m T ))a s {m T ) n (3.12) 



n=0 



f ds { s s 3 s A \ 1 
A n {a s {m T )) = J — I 1 - 2-^ + 2^ - J ^ —^—^ ln (_ s/m 2 )]n+ i > 

|s|=mj 



where d® are the expansion coefficients of Q 2 dU/dQ 2 , defined as in Eq. (|2.3|) . 
Mj} D (a s (m T )) and d®' MS are given in the last two columns of Table |IJ From this 
we would conclude from truncation of the series that M^ 13 is about 2.05 or even less 

since the coefficient <i^' MS = 3.10 is almost a factor three larger than the exact a 3 .- 
coefficient 1.26 for Q 2 dU/dQ 2 , which indicates that keeping vacuum polarization alone 
is not a good approximation^ for Q 2 dU/dQ 2 beyond order a 2 . 

The difference between M^ M5 (0.229) = 2.233 and 2.05 can be explained by the 
different weight of infrared and ultraviolet regions of internal momenta for R T as com- 
pared to Q 2 dU/dQ 2 . Since the A n in Eq. ( |3.12|) are positive numbers of order one, the 



behaviour of M^' LDP (a s (m T )) as function of N is controlled by the series for Q 2 dU/dQ 2 . 
The point is that the series for this quantity becomes dominated by ultraviolet regions 
of integration much earlier than R T itself, for which the contour integral suppresses the 
ultraviolet renormalon singularities as well as their residues. For instance, the ratio 
of the contribution of the leading infrared renormalon pole and the leading ultraviolet 
renormalon pole to the coefficients of the expansion for R T at order n + 1 is given by 

^•f(i)"n (3.13) 

and the crossover to the asymptotically expected ultraviolet renormalon dominance 
takes place at n w 8 — 9. For Q 2 dU/dQ 2 we obtain instead 

with crossover at n ~ 4 — 5, which is indeed confirmed by the coefficients for d^' MS 
given in Table |l|. The problem with the onset of ultraviolet renormalon dominance is 
that then one is forced to construct the analytic continuation of the Borel transform 

9 Thc discrepancy is partially removed, when two-loop running is incorporated, see Sect. 5.2. 
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in order to overcome the formal l/m^-uncertainty, associated with the truncation of 
the series due to ultraviolet renormalons. Since for a s = 0.32 this onset of divergence 
is expected at n ~ l/a s ~ 4 — 5, we see that, in the MS scheme, the resummation 
of Le Diberder and Pich can not be continued beyond n « 5 without running into 
this l/m^-uncertainty. It is the suppression of ultraviolet regions of integration that 
allows a fixed-order approximation of R T itself to higher n than five, which explains 

the proximity of Mj ,M (0.229) from truncation to 2.233 in this case. Recall that after 
resummation, there is no uncertainty left due to ultraviolet renormalons, because the 
sign-alternating behaviour is summable. 

Let us turn to the determination of a s . We may write the (normalized) hadronic r 
decay width as 

R r = 3 (\V ud \ 2 + \V US \ 2 ) S EW {l + + Sew + 5 P } , (3.15) 



where Sew an d Sew are electroweak corrections, S^°> is the perturbative QCD correction 
(with quark masses set to zero) and S p denotes power suppressed corrections in 1/m 2 ., 
including quark mass and condensate terms. We have borrowed the values for Sew, 
5ew and S p from |3T, [2^] . To a very good approximation we can neglect the variation 



of S p with a s and we have evaluated it at a s {m T ) = 0.32. Then the experimental value 
R T = 3.56 ± 0.03 (quoted from ||22||), obtained from an average of branching ratio and 
lifetime measurements translates into a constant experimental value^] 

6<®, = 0.183 ±0.010 (3.16) 



for the perturbative QCD corrections. In Fig. |3| we have plotted the prediction for 5^ 
as a function of a s (m T ) after resummation of one- loop running effects (recall for R T 
the exact a 2 - and a^-coefficients are very well approximated by one-loop running) as 
compared to the a 2 - and (^-approximation as well as the (^-approximation, including 
the resummation of ||10[ . We conclude that resummation of one-loop running reduces 
the central value for a s (m T ) by approximately 10% to 

a s {m T ) = 0.29. (3.17) 



There is some difficulty in assigning an error to this value, which is related to 
the extent to which one trusts the restriction to one-loop running effects as a good 
estimate of higher order coefficients, not known exactly at present. In view of the 
fact that one-loop running overestimates the coefficient d 2 of Q 2 dU/dQ 2 , one might 
consider M^ AS (a s ) as an upper estimate for the effect of higher order perturbative 
corrections and the quoted value for a s {m T ) as a lower estimate for the central value. 
On the other hand, we shall see in Sect. 5.2, that a partial inclusion of two-loop running 
points towards even a lower value of a s (m T ). It seems safe to conclude that higher 
order corrections add up constructively and a cumulative effect is most likely to shift 

10 The latest experimental data, which we have not taken into account, seem to indicate a larger value 
for R T . In this case, too, resummation reduces a s (m T ) obtained without resummation by 10% — 15% . 
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°' 1 0.26 0.28 0.3 0.32 0.34 0.36 0.38 



Figure 3: Perturbative corrections to the r hadronic width. I: After resummation 
of one-loop running effects; II: Exact order (^-approximation; III: Exact order a^- 
approximation including the resummation of running coupling effects along a circle in 
the complex plane (taken from Pich); IV: Exact order (^-approximation for comparison. 
The shaded bar gives the experimental value with experimental errors only. 



a s (m T ) towards the above value. Optimistically, one could even hope for a reduction 
of the theoretical uncertainty from the unknown residual perturbative corrections. A 
conservative evaluation would not exclude the entire region from 0.27 to 0.34 for a s 
at the scale m T . The lower bound reflects experimental and other than perturbative 
theoretical uncertainties and the upper bound follows from an analysis that takes into 
account only the completely known perturbative corrections up to cubic order. 

The implications of the expected ultraviolet renormalon dominance for R T in large 
orders, which we have briefly touched upon above, have recently been investigated in 



Ref. p7| ) whose authors use conformal mapping techniques to construct the analytic 
continuation of the Borel transform beyond its radius of convergence set by the first 
ultraviolet renormalon. As noted in [2jJ, such a technique is not particularly successful, 
when the series is not already close to the asymptotic ultraviolet renormalon behaviour 
(which is not the case for R T and Q 2 dU/dQ 2 up to cubic order), because then the 
intricate cancellations necessary to push the ultraviolet renormalon further away from 
the origin of the (conformally mapped) Borel plane than the first IR renormalon do not 
take place. The variation of results obtained from different mapping functions can then 
be considered analogous to the variations induced by different choices of renormalization 
schemes and as in the latter case it is difficult to decide to what extent such a (in 
principle arbitrary) variation should be considered as a theoretical error. Based on 
the evidence presented above, we believe that ultraviolet renormalons are sufficiently 
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suppressed in fixed-order perturbative approximations to R T (but potentially not to 
Q 2 dU/dQ 2 ) to be safely ignored at present. 



3.2 (AQCD/^r) 2 -corrections 

In this Subsection we investigate whether resummation of perturbative corrections can 
introduce power corrections, which elude the operator product expansion. We will ex- 
clude from the discussion the effect of renormalons, which have received main attention 
in this context p5| , |TT| , |2"P| , |T3| |5U| , As far as evidence from explicit calcula- 



tions is available, infrared renormalons are in correspondence with condensates in the 
OPE as expected. There is no indication for explicit power corrections of dimension 
two from this source, which could turn out to be numerically significant. The effect 
of the dominant ultraviolet renormalon divergence is taken into account automatically 
by the error estimate due to unknown higher order perturbative corrections. After 
resummation it disappears completely in principle]^], although in practice this is not 
simple to implement [^J without approximations (like large fio). It is also conceptually 
important to distinguish the statement that 1 / Q 2 -corrections are absent in the OPE of 
correlation functions at euclidian momenta from that that 1/m^-corrections are absent 
in R T . Validity of the first needs not imply the second, though the second will hardly 
be true without the first. 

3.2.1 Definition of the coupling 

Our first point of concern is the definition of the coupling parameter inherent to pertur- 
bative expansions and therefore to their resummations. To emphasize that this question 
is not connected with the minkowskian nature of R T , we work with the derivative of 
the correlation function D(Q 2 ) rather than with R T itself. 

Let us first take a closer look at the Borel sum (to be definite, let us assume a 
principal value prescription throughout this subsection, which corresponds to taking 
the real part of r (X 2 L )) in the representation as an integral over finite gluon mass in 
the large-/?o (NNA) approximation. Let us denote the two terms on the right hand 
side of Eq. ( p.24|) by I(a s ) for the integral and L(a s ) for the contribution from the 
Landau pole in the dispersion relation for the running coupling (using the technique 
of Sect. 2.3 or Q). The leading (two-loop) radiative correction with finite A can be 
obtained from the Borel transform by a Mellin transformation, see Eq. ( |2.16| ). Taking 
the integral analytically can be difficult, but for our present purpose we are interested 
only in the coefficients in the small-A 2 expansion. This expansion can be obtained 
easily, since particular terms proportional to (X 2 /m 2 ) n (modulo logarithms) correspond 
to contributions from singularities of the integrand in the right li-plane at u — n. For 
example, to pick up the term of order X 2 /m 2 one can evaluate the integrand in Eq. ( 2.16 ) 
at u — 1, except for T(—u). For the D-function defined in Eq. (|3.8| ) we obtain 



11 This is seen explicitly in the representation Eq. ( 2.24 ) for the Borel integral, which avoids con- 
struction of the analytic continuation of the Borel transform beyond its radius of convergence set by 
the first ultraviolet renormalon. 
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+ 0(A 6 ln 2 A 2 /g 2 



JC(3) 



A 2 
Q 2 



21n(Q 2 /A 2 



20 
T 



K(3) 



A*_ 
<3 4 

(3.18) 



The presence of quadratic terms, X 2 /Q 2 , is not in conflict with the operator product 
expansion, since such terms come from the region of large momenta: As emphasized 
in Sect. 2.4, only non-analytic terms in the small A 2 -expansion can unambiguously be 
identified with infared contributions. The leading non-analytic term is proportional 
to A 4 In A 2 and produces a correction proportional to l/Q 4 , which can be related to 
the contribution of the gluon condensate |3l| in the OPE. This contribution agrees 
with the calculation of the gluon condensate with finite gluon mass in p2[ , after the 
corresponding Wilson coefficient is extractedf 2 } (The difference in the constant 20/3 
arises, because we consider the derivative of II {Q 2 ).). 

The presence of a A 2 /Q 2 -correction in Eq. ( |3.18|) implies that the Landau-pole con- 
tribution in Eq. Q2T211 ) is 



L(a s (Q)) oc exp 



oc 



A 2 

QCD 

Q 2 



(3.19) 



Numerically, this term is quite substantial. Taking a s (Q) = 0.32 (which corresponds 
to Q ~ m T ), separation of the two terms in Eq. (|2.24j ) for D(Q 2 ) amounts to 



M^(0.229) = 1.48 = integral over gluon mass + 0.31 . 



(3.20) 



Note that L(a s ) has identically vanishing perturbative expansion. Thus, without any 
additional information, keeping I(a s ) alone in Eq. ( |2.24j ) provides an equally legitimate 
summation of the original series, which differs from Borel summation by terms of order 
1/Q 2 , which are not related to renormalons or any particular regime of small or large 
momenta. We conclude, at this stage, that statements about power corrections are 
meaningful only with respect to particular summation prescriptions. 

Physically, dropping the contribution L(a s ) is equivalent to dropping the Landau 
pole contribution to the dispersion relation for the running coupling which can be 
interpreted as a redefinition of the coupling, such that the new coupling has no Landau 
pole. Since Borel summation in our limit of large (3q coincides (in the V-scheme, C=0) 

12 The full contribution of order A 4 /Q 4 which is multiplied by ln(Q 2 /A 2 ) + const should be schemat- 
ically decomposed as ln(/z 2 /A 2 ), contributing to the gluon condensate, and ln(Q 2 /^ 2 ) + const, con- 
tributing to the coefficient function to a s accuracy. Here \i is the scale separating small and large 
distances. If one considers A as as an infrared cutoff, which is natural, then the full contribution 
should be ascribed to the coefficient function. Hence the rationale of ascribing the constant term to 
contributions of large momenta (small distances). In general, the separation of matrix elements and 
coefficient functions is of course factorization scheme-dependent and constant terms can be reshuffled. 
We see once more that only non-analytic terms (logarithmic in this example) can be used to trace 
infrared contributions. 
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with averaging the running coupling a s (k) = a s (Q)/(l— /3 a s (Q) \n(k 2 /Q 2 )), it is readily 
seen from Eq. (|2.33|) that neglecting L(a s ) corresponds to averaging with a coupling 
af £ (Q), related to ot s {Q) by 

af (Q) = a s (Q) - - ^ = a s (Q) - 1 e V(/WQ)) + ( e V(/WQ))) . ( 3 .21) 

This coupling has no Landau pole and freezes to a finite value as Q 2 approaches zero. 
However, it has l/Q 2 -corrections to its evolution and correspondingly to the /^-function: 

eS (af) = Maf) 2 + 2af ) e VCA>af ) + ^/(ftaf)) (3 22) 

Note that the absence of a Landau pole in off is seen only after summation of an 
infinite number of power corrections in 1/Q 2 (exponentially small terms in a s (Q)) in 
Eq. ( |3.21| ). Further, the coefficients of perturbative expansions of any quantity are the 
same whether one uses a s (Q) or af(Q) and diverge although the average 

J d 4 kF(k,Q)^^- ~ I(a s (Q)) (3.23) 

has no Landau pole ambiguities. We emphasize once more that averaging one-loop ra- 
diative corrections with this freezing coupling differs from Borel summation of the series 
in a s (Q) (which gives I(a s (Q)) + L(a s (Q)) by l/Q 2 -terms. So does Borel summation 
of the identical series in a1 s (Q), because the couplings differ by such terms. We shall 
argue that couplings like off(Q) obscure the relation to the operator product expansion 
in the sense that explicit l/<5 2 -terms must be added to the resummed result (such as 
L(a s (Q)), had one used a^ ff (Q)) in order to cancel spurious l/Q 2 -effects in the large 
Q 2 -expansion of the resummed result. This remark applies identically to a freezing 
coupling of type l/(— Ajln(c + Q 2 / A 2 )) and thus to the procedure of 0, where such 
a coupling has been used to estimate the size of infrared contributions for quantities 
related to heavy quark expansions. 

We will now try to make precise the statement that l/Q 2 -terms should be absent in 
the OPE. In order to talk about power corrections we have to attach some meaning to 
the divergent perturbative expansion at leading order in 1/Q 2 . A natural (though one 
might still ask whether it is justified to prefer Borel summation- type schemes to any 
other) definition apparently is 

£(Q 2 ) ex a Ct - BS[D(a s (Q)) pcrt }\ - O (Jp) , (3.24) 

where BS[D(a s (Q)) peT t\ denotes the Borel integral (with principal value prescription) 
of the perturbative expansion of D in some a s (Q). However, this statement is still 
ambiguous, because it implies knowledge of the Q 2 -dependence in the coupling, which 
is arbitrary to a large extent. Moreover, it is not sufficient to appeal to the usual 
ambiguities in the choice of perturbative renormalization schemes, since the coupling 



29 



and its evolution must be specified to power-like accuracy. We can bypass this point, 
noting that QCD with massless fermions has only a single free parameter. Since we 
are discussing asymptotic expansions in a dimensionful parameter 1/Q 2 , it seems most 
natural to choose a physical mass scale as this parameter, say m 2 . This is especially 
natural in the context of lattice definitions of QCD, where one would trade the bare 
coupling for m 2 . Then, if the operator product expansion exists nonperturbatively, this 
suggests the existence of a double expansion in 1 / ln(Q 2 / m 2 ) and m 2 /Q 2 at large Q 2 
and the interpretation of the statement that there are no 1/Q 2 terms could be0 

(TYl 2 \ ^ 
' (3 " 25) 

Here BS[D(Q 2 ) peit ] denotes the Borel sum (with principal value prescription) of the 
leading term in the expansion of D(Q 2 ) at large Q 2 , which is an infinite series in 
1/ hi(cQ 2 /m 2 ), where c is a constant to be specified later. The constant K(8, R) depends 
on the opening angle 5 and radius R of the sector in the complex l/Q 2 -plane, where 
Eq. (|3.25|) is supposed to hold. In general, K(5, R) will neither be continuous nor 
bounded as a function of these two parameters. In particular, one can not expect a 
uniform bound in the entire cut plane. In the limit 5 — > n, this is related to violations 
of duality. We stress that as far as mathematical rigour is concerned, the validity 
of Eq. ( p.25| ) must be regarded as purely hypothetical. We wish to present it as a 
mathematical formulation of the assumption that the operator product expansion holds 
at euclidian momenta (i.e. 5 > 0) and no l/<5 2 -terms are present in the asymptotic 
expansion. We note that the condition Eq. ( |3.25|) is stronger than the condition that 
long-distance contributions can be factorized into condensates, which does not exclude 
the presence of power-like corrections, in particular 1/Q 2 , to coefficient functions, in 
particular of the unit operator, from short distances. 

After we have chosen m p as the fundamental parameter of QCD, we may define 
the coupling by its beta-function and an overall scale. We define (again, to leading-/? 
accuracy) 

asiQ) = -PoHcQi/mj) = -PoHQ 2 /A 2 v ) ' (3 - 26) 

where we have fixed c by matching the large Q 2 behaviour with the V^schemeQ It is 
easy to see that Eq. ( |3.25| ) implies Eq. ( |3.24j ) with this coupling, which by definition 

13 This is still a simplification. In fact, lnln(Q 2 / TO p) is also expected to appear. We will restrict 
the discussion to the large-/3o limit, which might be considered as an analytic continuation to a large 
negative number of massless fermion flavours. In this limit, only ln(Q 2 / 'm 2 p ) appears and Borel trans- 
formation with respect to ln(cQ 2 / 'm 2 p ) has a unique meaning. 

14 This is a matter of convenience, since it eliminates writing C in the large-/3o approximation. We 
could also have matched to the perturbative MS coupling. We also note that the change of variables 
from m„ t o a s {Q) is singular due to the Landau pole at A v . However, this is not a restriction, since 



Eq. (3.25) is limited to finite R anyway. The position of the pole may be varied by the choice of c, 
which implies a reorganization of powers in inverse logarithms, but the pole occurs always at a finite 
value of Q 2 . 
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has no power-like evolution. This implication is not valid for couplings, that incorpo- 
rate l/Q 2 -dependence in their running. It is in this sense that we believe the use of 
freezing couplings is hazardous. //Eq. ( |3.25| ) is correct, then the use of Borel summa- 
tion for perturbative series expressed in terms of such a coupling, or averaging lowest 
order radiative corrections with such a running coupling, necessitates the addition of 
explicit l/Q 2 -corrections simply to cancel such corrections hidden in the definition of 
the coupling. 

In practice, in one way or another, one relates physical quantities and an unphysical 
coupling like Eq. ( |3.26p might be considered as an intermediate concept only. How- 
ever, the importance of being definite with the evolution of the coupling to power-like 
accuracy is not diminished by the use of physical couplings. To give a somewhat con- 
structed example: If one expressed R e + e - as an expansion in the effective coupling, 
defined by QCD corrections to the Gross-Llewellyn-Smith sum rule, one would expect 
l/Q 2 -corrections to this perturbative relation, which are imported from the definition 
of the coupling. 



3.2.2 Analyticity and the Landau pole 

After renormalization group improvement, perturbative expansions are plagued by the 
Landau pole. This unphysical singularity is endemic not only to perturbative expan- 
sions, but to the operator product expansion, truncated to any finite order. It requires 
some care in defining resummations for quantities like R T , which are related to euclidian 
quantities by analyticity. For the purpose of illustration, we shall restrict ourselves to 
the approximation, where P(a s ) = Poa 2 and adopt the V-scheme, C = 0, in all explicit 
formulas. Then, ignoring an irrelevant Q 2 -independent subtraction, we can write the 
Borel sum (with principal value prescription, as usual) as 



BS[U](Q 2 ) = J due - u ^^B[U](u,Q 2 / f , 2 ) = ^Jdui^\ F(u) 

^ ' 

= BS[U] < (Q 2 ,u ) + BS[n} > (Q 2 ,u ), (3.27) 



where we use a s (Q) = l/(— fto ln.(Q 2 /A v )) and F(u) is Q ^independent. In the second 
line, we have defined two new functions by splitting the Borel integral into two regions 
from to u Q and u to 00. From the explicit form in Eq. (|3.9Q , we deduce that, for any 
(positive) Uo, BS[H] < (Q 2 , u ) is analytic in the Q 2 -plane cut along the negative axis. 
On the other hand the it-integral in BS[U}(Q 2 ) diverges, when Q 2 < Ay and Q 2 = Ay 
is a singular point. We say that BS[U](Q 2 ) has a Landau pole, though, in general, 
Q 2 = Ay will rather be a branch point. 

Notice, since in practice one truncates the operator product expansion at operators 
of some dimension d, it is perfectly consistent to replace BS[U](Q 2 ) by BS[Il] < (Q 2 , uo), 
provided we choose uq > d/2 + 1, since the difference is bounded by 



BS[U](Q 2 )- BSiUUiQ 2 ,^) <K(R)(-£) (3.28) 
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and vanishes faster than other corrections neglected in the truncation of the OPE. In this 
case, the bound can be established in a cut circle of radius R in the l/Q 2 -plane. Again, 
it will not be possible to establish uniform bounds, since typically K(R) ~ 1/ \n(RAy). 
Thus, although for fixed Q 2 > 1/R, the difference can be made vanish faster than any 
desired power of Ay/Q 2 by increasing Uq, at fixed uq the bound may become arbitrarily 
weak as Q 2 approaches Ay. Still, we conclude that the presence or absence of a Landau 
pole in resummed results is related to the behaviour of the Borel integral at infinity 
and is thus an effect that formally vanishes faster than any power of 1/Q 2 - 

Consider now the two different representations for the tau decay width, Eqs. (|3.3|) 
and (|3.5| ), in this light. We use equality of vector and axial- vector correlators in per- 
turbation theory and abbreviate the weight function by w(s/m 2 ). Then 

K e q .(3. 3) = 24?r f *l w (j_\ AdiscssrnKs), (3.29) 

J to- \miJ 2i 




(3.30) 



Note that in both cases summation is carried out inside the s-integral. The following 
considerations can easily be extended to the situation, where summation is taken after 
s-int£gration (and donot lead to any of the differences observed below). We define 
R^f3(u ) and R^f^(u ) by the replacement of BS[U](s) by BS[n]<(s,« ). Using 
the analyticity properties discussed above as well as m 2 > Ay, it is straightforward to 
find 




(3.31) 



but 

^.(3.3) _ ^,0) = _ U7ri f *L W (*) BSiUUs,^) , (3.32) 

J mi. \mi-J 
c 

where the contour C runs along a circle of radius A v around s = —Ay. Against 
appearances, the right hand side is independent of u . Of course, one can not take u$ 
to infinity inside the integral and conclude that it is zero. BS\U\ > {s ) uq) has a pole or 
branch point at s = —A v and we find 

^.(3.3) _ ^eq.Q _ A| 

and similarly 




(3.34) 
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Since Eq. ( |3.28| ) is not valid for s < A v , this result should not surprise. The difference 
in Eq. ( p.33|) arises, because the resummation introduces (or preserves) the Landau 
pole singularity, which is in conflict with the analytic properties of the exact correlation 
function, that have been assumed to derive Eq. (|3.5|) from Eq. (|3.3|) . 

Should one conclude then that resummations of perturbative corrections are ambigu- 
ous by terms of order even if there are no 1/Q 2 terms in the OPE at euclidian 
momenta in the strong sense of Eq. ( |3.25| )? Do we have evidence for power corrections 
not captured by the OPE, since they originate from u = oo in the Borel integral? The 
answer is no. A positive answer would be warranted, if there were no reason to prefer 
the prescription Eq. ( |3.29|) to Eq. (|3.30|) or vice versa, while only the second can be 
used. The reason is that the region \s\ < A v can not be penetrated to any finite order 
in the short distance expansion (l/Q 2 -expansion). Since all summation prescriptions 
of perturbative expansions are formulated within the context of the short- distance ex- 
pansion^], they can not be applied to \s\ < A v , unless the summation of this expansion 
itself is understood. This discards Eq. ( 3.29 ) as a legitimate summation. One must first 
use the analyticity properties of the exact correlation functions to deform the contour 
outside the region \s\ < A v , before an attempt at summing perturbative expansions can 
be made, which privileges Eq. (|3.5| ) as the starting point. A different way to express 
this fact is to observe that the principal value Borel integral is defined for \s\ < A v only 
in the sense of analytic continuation, but cannot be used as a numerical approximation 
since all power corrections in the operator product expansion are of the same order 
of magnitude in this region. It is only when all these are taken into account that the 
Landau pole vanishes in physical observables. 

To conclude this Section, let us mention that a relation between analyticity in s and 
behaviour of the Borel integral at infinity has been noted in a very different and much 



more physical context in |j33| . The presence of resonances and multiple thresholds on the 
physical axis was observed to be in conflict with Borel summation^. In this case, the 
restoration of the correct analytic properties should also be understood in connection 
with summing the OPE itself p5| . A much more elaborate argument has been presented 
in P6 |, where the presence of resonances was connected with the divergence of OPE, 
which also offers a way to understand the concept of duality and its limitations. If the 
OPE is by itself only an asymptotic expansion, it is quite possible that its application 
is limited to a finite phase range in the complex s-plane around the negative s-axis. We 
believe that this possibility should be taken seriously, since it might imply the presence 
of 1/so-corrections in finite averages up to so of the discontinuity along the physical 
axis (for R T , sq = m 2 ). Unfortunately, we do not know how to substantiate or disprove 
such a statement theoretically. 



15 See the discussion below Eq. (p.S|). 

16 Historically, it is interesting to note that the horn-shaped analyticity region in the coupling, which 



leads to this conclusion, was discovered for the photon propagator in Q . 
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4 The pole mass of a heavy quark 



Although on-shell quarks do not exist and there is presumably no natural nonpertur- 
bative definition of a pole mass of a quark, the pole mass has proven useful as an 
auxiliary concept in applications of perturbative QCD to heavy quark physics, where 
physically the quark is expected to be close to its would-be mass-shell. Still, the relation 
between the pole mass and an off-shell renormalized mass is known to have large per- 
turbative coefficients from small loop momenta, at least in high orders of perturbation 
theory Jl7|, [57). Since neither quark mass definition is physical, this might appear as 



an irrelevant problem. However, the behaviour of perturbative expansions of quantities 
involving quark masses change accordingly with the quark mass definition used and 
in general one expects coefficients to be significantly reduced, when the pole mass is 
abandoned in favour of an off-shell renormalized mass |37], |2(J, such as MS. In the latter 



case, it is quite well-known that the exact two-loop coefficient [38| in the relation to the 
pole mass is substantial and one might wonder whether this is coincidental. 

In this Section we apply resummation and NNA to the difference between the pole 
and MS mass. Apart from its practical interest, we also use this quantity to illustrate 
certain features of higher order perturbative corrections, when the masses of fermions 
in loops are finite. The reader interested only in results may jump directly to Sect. 4.3, 
where our best estimate is presented. In this Section always refers to the MS 

coupling. 

4.1 Preliminaries 



To begin with, let us quote the exact two- loop result from [35]: 



5m 
m 



m 



pole 



m Ms( m Ms) = 4 a s( m Ms) 
3 7T 



m Ms( m Ms) 



1 + 4.68 {-(3 



N, 



0.89 
0.91 



a s( m MS") 



(4.1) 



The upper entry in curly brackets refers to the calculation with Nf massless quarks and 
a quark of mass m in loops, the lower entry to the situation, where the quark flavour 
of mass m is excluded from loopsQ. The superscript on (3 indicates the value of Nf 
to be taken for the /3-function. Since (-Po f ) ~ 2/3 for cases of interest, Nf = 3,4, 

Nt 

we note that keeping the term proportional to f3 alone provides indeed a reasonable 
approximation within 30 — 40% of the exact coefficient. 

Note that above we have normalized at the scale since we prefer to 

eliminate m po i e in all places. When Eq. (|4.1| ) is expressed in terms of the MS mass 
normalized at m po i e , 



^polc 



m- 



< M g(m p ole 



m 



MS 



(mpoie) 



4 a s (m P oic 
3 



7T 



1 



4.68 (-f3 



Nf 




a s {m 



pole J 



(4.2) 



17 To obtain the lower value, it is necessary to insert a missing Cf in front of the mass correction in 
Eq. (17) of @, see also @. 
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the approximation of keeping only the term proportional to (3 is significantly im- 
proved. This ambiguity is in fact a source of trouble within the BLM and our extended 
prescription. For Sm additional renormalization scheme and scale dependence is present 
from the definition of the quark mass parameter, whereas the BLM method by construc- 
tion deals only with the scale-dependence of the coupling. Although the calculation of 
fermion loop insertions and subsequent restoration of /3q also provides the anomalous 
dimension of the quark mass within the same approximation, the scheme ambiguity in 
the size of neglected genuine two-loop corrections is amplified by this additional source 
of scheme- dependence. If ji\ — ji 2 . ~ 0(a s ), then the difference mj^(fii) — mj^g(/i 2 ) 
must be neglected in the approximation of large /3o, although it can be sizeable, if the 
coefficient of a s in n\ — /j% is so. In the following, we will work with Eq. ( }4.1|) , whenever 
relevant. 

Let us first consider the (excellent) approximation, that the quark flavour of mass 
m is neglected in loops. All other Nf quarks are massless. In this case the relevant 

Nt 

/3-function coefficient is f3 and an exact expression for the Borel transform of the mass 
shift exists |17||: 



B [5m/m] (u) = — 
3ir 



5«/3 P M ..N r («) r (l-2U) , G {U) 



e^6(l-n) ^ r(3 - ' + 



(4.3) 



Gq(u) is defined in the following way (see Appendix A): If g n are the expansion co- 
efficients of Gq(u) in u, then Gq{u) has expansion coefficients g n /n\. Gq{u) can be 
calculated with the methods of Appendix A and is found to be 

GoW = 4(3 + 2, ) f(r -^j|__ . (4 . 4) 

Higher order perturbative corrections with fermion loop insertions into the one-loop 
diagram can then be computed according to Eq. (|2.7| ). Equivalently, one can start from 
Eq. ( |2.28| ) with the input of the one-loop mass shift with finite gluon mass, given by 
(x = \ 2 /m 2 ) 



ro(A 2 
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3^ 



4 + x In x 
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\/4 - x 



arctan 
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l + ^ + l x + (x^ 2 ) 



(4.5) 



We define coefficients d n as in Sect. 2 by 

5m AaJm 



m 3 7r 



n=l 



(4.6) 



They are listed in Table 0. Higher order perturbative corrections grow very rapidly, 
as expected from the dominant asymptotic behaviour from the pole at u — 1/2 in 
Eq. 
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n 




M n [c- quark] 


M n [b-quark] 





1 


1 


1 


1 


4.6861511 


2.176 


1.623 


2 


17.622650 


3.286 


1.935 


3 


109.85885 


5.024 


2.193 


4 


873.92393 


8.492 


2.467 


5 


8839.6860 


17.30 


2.835 


6 


105814.28 


43.76 


3.420 


7 


1484968.4 


137.0 


4.514 


8 


237407365 


511.0 


6.838 


oo 




1.712 ±0.608 


2.041 ±0.201 


* 




0.096 m c 


0.096 m b 






0.437 m c 


O.U7m b 



Table 2: Higher order perturbative corrections to the difference between the pole and 
MS mass from n fermion loop insertions. The second and third columns give partial 
sums and summations for the charm quark [a s = —j3 a s (m c ) = 0.251) and bottom 
quark (a s = — (3Q^a s (m b ) = 0.133). In the latter case, the c-quark is treated as massless 
inside loops. 



d n e 5 / 6 2 n n\ . (4.7) 

This asymptotic formula is in fact a very good approximation (within 5%) to the exact 
coefficient d n for all n > 1, which seems to imply that the saddle point approximation 
of the loop momentum distribution inherent to deriving the asymptotic behaviour is 
a very good substitute for the exact distribution, even if the width of the Gaussian 
is not small. Given that genuine two-loop corrections are rather small compared to 
di, it can be asserted that the exact two-loop coefficient is already dominated by the 
first infrared renormalon at u = 1/2. Evidently, this observation is scheme-dependent 
and we do not know why the MS scheme is preferred. Let us note, however, that a 
similar coincidence can not be expected and indeed does not happen (see Sect. 3) for 
quantities, whose asymptotic behaviour is dominated by ultraviolet renormalons, like 
R e + e - or R T , since with C = —5/3 in the MS scheme, the leading ultraviolet singularity 
at u = — 1 is suppressed compared to the leading infrared singularity at u = 2 by a 
factor e -5 ~ 7 ■ 10 -3 . Therefore the latter is expected to dominate in intermediate 
orders JJIJ (if there is any regularity at all) and it might not be an accident, that exact 
low order coefficients are indeed of same sign for these quantities. 

Table also shows how the one-loop radiative correction to the pole mass is mod- 
ified by inclusion of a finite number of fermion loops and by summation according to 
Eq. ( |2.29p . We have taken a s {m b ) = 0.2 and a s {m c ) = 0.35. Recall that the sum 
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corresponds to a principal value prescription for the Borel integral. The errors quoted 
correspond to the imaginary part of the Borel integral, when it is defined by deforming 
the contour into the complex plane. The imaginary part can be taken as an estimate of 
inherent uncertainty of perturbative relations. We have actually divided this imaginary 
part by n, which upon inspection we find closer to, but somewhat smaller than the 
naive estimate of uncertainty by the minimal term of the series. A more conservative 
estimate would be to enlarge these errors by a factor of two. The BLM scales in Table ^ 
are defined by (cf. Eq. (|1.5| )) 



The uncertainty in Moo(a s ) translates into an uncertainty of the BLM scale m*^ by the 
previous equation. We want to point out that the usual BLM scale m*, which uses 
only di as input is smaller than although all higher orders d n add up positively. 
The reason is that a s (ra*) upon re-expanding in terms of a s (m) implies d® LM = d™. 
For the charm and bottom quark mass the most important effect comes from n = 1,2. 
But since d± is rather large, <if LM = d\ « 21.9 > 17.6 = d 2 , and the usual BLM 
prescription overestimates the size of radiative corrections. This behaviour is quite 
generic to quantities dominated by scales below a few GeV and with a leading infrared 
renormalon at u = 1/2. 

4.2 Effect of internal quark masses 

Before proceeding to realistic charm and bottom quarks we want to illustrate the effect 
of finite quark masses inside loops on the coefficients of diagrams with quark loop 
insertions in the toy example of the mass shift due to a single quark flavour with mass 
rrii in loops. The "heavy" quark of mass m is again excluded from loops. 

The factorially large contribution from small and large loop momenta can be traced 
to the logarithmic behaviour of the vacuum polarization for a massless particle at very 
small and very large virtuality. In the case of small virtuality, the large coefficient arises 
from momenta k ~ me~ n <C m. For a massive quark the logarithmic behaviour is cut 
off, because at very small momenta its vacuum polarization is proportional to k 2 /mf. 
Thus, no matter how small the quark mass, the factorially large contribution should be 
eliminated when n is such that me~ n < rrii. There are no infrared renormalon singu- 
larities in the Borel transform in the absence of massless particles. This expectation is 
verified in Fig. |], where we have plotted the ratio of the coefficient d n (rrii), computed 
with internal mass rrii, and d n (0) for the massless case as in Table |2] as a function of the 
number of loops n for various ratios of rrii/m. Since rrii serves effectively as an infrared 
cutoff, Fig. |] provides some information on what proportion of the coefficient d n (0) 
originates from low momentum regions. Note that low momentum means low momen- 
tum compared to m and not Aq CD . Eventually, as n becomes large, coefficients will be 
dominated by large momentum for any rrii and the asymptotic behaviour of Eq. (|4.7| ) 
is replaced by a sign-alternating behaviour due the leading ultraviolet renormalon 



ml = m exp -- — (M\(a 8 )-l) 



s 




(4.8) 
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Figure 4: Ratio d n (rrii)/d n (0) for different values of internal quark masses rrij as a 
function of the number of fermion loop insertions. The value of m 2 /m 2 is indicated to 
the right of each curve. 



d n (m.i 



n>l 



-1 + 



9 m\ 
2^2 



-5/3 



(4.9) 



We find that for m 2 /m 2 > 0.1, this asymptotic behaviour practically coincides with the 
exact coefficient for n > 6. 

The presence of mass dependence in the coefficient of ultraviolet renormalons re- 
quires some explanation. The singularity of the Borel transform at u = —1, which is 
responsible for Eq. Q4.9| ), is due to the presence of a l//c 6 -term in the expansion of the 
Feynman integrand for the one-loop mass-shift at large k: 



5m ~ a, m I d k 



a m 

¥ + ~k? 



+ ... 



(4.10) 



The l//c 4 -term creates a logarithmic ultraviolet divergence, which is subtracted by min- 
imal subtraction. The effect of the Borel transformed gluon propagator (with fermion 
loop insertions) amounts to insertion of (see Sect. 2.3) 



exp(— vXl(k 2 ) / a s ) 



m 



-c 



Tfl- 

l + 6u^ + 




(4.11) 



into the integrand. The pole at u — — 1 arises, because close to u = —1 the l/fc 6 -term 
is converted into a logarithmically divergent term. But when rrii is not zero a second 
contribution of order 1/k 6 arises, when a/fc 4 combines with the first mass correction 
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in Eq. ( 4.11 ). This term accounts for the mass-dependence of ultraviolet renormalons 
in minimal subtraction schemes. Had one chosen a renormalization prescription, which 
subtracts the a//c 4 -term inside the integrand, that mass-dependence would be absent for 
the first ultraviolet renormalon at u = —1, but still present in coefficients of singularities 
at u — —2, — 3 . . .. 



4.3 Charm and bottom pole masses 

For realistic charm and bottom quarks, the numerical results of Sect. 4.1 need to be 
amended in two respects: The quark, whose mass shift is considered, should also be 
taken into account in loops. This effect is tiny, but we include it for completeness. 
More interesting is the effect of finite charm mass on the bottom pole mass, because 
the charm mass might be larger than the typical loop momentum already in low orders. 
If so, it would seem more appropriate to use (3^ rather than [3^, when restoring the 
QCD /3-function coefficient from the calculation of fermion loops. Strictly speaking, 
changing the mass of one flavour is a negligible effect in the formal large-/3o limit (i.e. 
large- Nf), and should be discarded in a consistent approximation. In realiter the effect 
is numerically noticeable and we prefer to include it as a calculable correction, keeping 
in mind its size as an uncertainty in our approach. 

To gain some understanding of the numerical results to follow, we trace the decou- 
pling of internal charm loops for the bottom pole mass in more detail. For this purpose, 
we ignore again the bottom quark in loops. When the order of perturbation theory is 
sufficiently large, internal integrations are dominated by momenta much smaller than 
m c . Let us denote coefficients including charm loops by d^ +1 \m c ) (— fi^)'"' and those 
without charm loops by d$ (—/3^) n . Then one might expect 

4 3+1 l(m c ) (-4 4) ) n a s (m b ) n+1 « cf a s (m h ) n+l (4.12) 



if n is sufficiently large up to corrections suppressed by a power of the typical internal 
momentum divided by m c . This is not quite correct, because there is no manifest 
decoupling in the MS scheme. In the limit that we consider only the contribution from 
small A 2 , appropriate for sufficiently large n, the Borel transform including the massive 
charm in Eq. ( [2.371 ) reduces to 



(4.13) 



up to power corrections. If the first factor were absent, this Borel transform would 
coincide with the one for purely massless flavours up to a rescaling of u which is just 
what is necessary to obtain Eq. ( |4.12[ ) from Eq. (|2.7| ). The first factor is present, because 
the vacuum polarization for charm has been renormalized in the MS scheme and not by 
zero momentum subtraction, which would lead to manifest decoupling. The difference 
can be accomodated by a change of the coupling constant. This is most easily seen 
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n 


4 3] 


43+1] 


4? (flSW) B 


M[ 3+1 l 

n 





i 


1 


i 


1 


1 


4.6862 


5.0984 


5.0613 


2.183 


2 


17.623 


20.514 


20.555 


3.291 


3 


109.86 


138.58 


138.38 


5.021 


4 


873.92 


1188.3 


1189.0 


8.470 


oo 








1.720 ±0.608 


m* 








0.078 m c 










0.406 m c 



Table 3: Perturbative corrections to the pole charm quark mass with (column 3) and 
without (column 2) charm inside loops. For comparison, column 4 gives coefficients in 
the limit that all internal masses are large compared to the typical loop momentum. 
The last column updates the modification of one-loop corrections due to coefficients 

4 3+1] - 



by combining the extra exponential factor with the exponential in the Borel integral, 
Eq. Q2.8D . We deduce that (up to power corrections in the typical loop momentum 
divided by m c ) 

(-P^) n a^(m b ) n+1 « 4 31 (-4 3 T «i 3) M n+1 , (4-14) 

where the two couplings are related by 

1 1 l_ j£_ 

ai 3 V)~ ai 4 V) + 67r 

Alternatively, we can simply replace u in the extra exponential factor in Eq. (|4 . X 3|) by 
the location of the closest infrared pole of the Borel transform, u = 1/2, and obtain 

dl +l \m c ) (-$ ) ) n a s (m b ) n+1 « exp (—L^\n^j d® a s (m b ) n+1 , (4.16) 



which is equivalent to Eq. ( 4.14| ) for large 



n. 



Numerical values for corrections to the charm pole mass and the bottom pole mass, 
including the charm and bottom masses in loops, are given in Tables |3| and [|. For the 
bottom pole mass, the result depends only on the ratio m c /m b and we used m 2 c /ml = 
0.1. For comparison coefficients according to Eq. (|4.16j ) are given. The last column in 
both tables displays the modification of the one-loop mass shift by running coupling 
effects. Note that at the charm scale, the perturbative series has to be truncated already 
at second order. The lower rows give the result of summation according to Eq. (|2.37Q 
and the corresponding BLM scales, Eq. (|4.8|) . Two technical observations are in order: 
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4.6862 


5.3093 


4.9774 


1.648 
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17.623 
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1.986 
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109.86 


159.69 


160.82 


2.276 
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873.92 


1.5 • 10 3 


1502.0 


2.601 


oo 








2.099 ±0.224 


m* 








0.064 m b 










0.117m& 



Table 4: Coefficients as in the previous table, where rf[f +2 ' now includes charm and 
bottom masses with m 2 /m 2 = 0.1. K = exp(l/(l2ir[3^) ln(m 2 /m 2 )) ps 0.9. 



First, as mentioned in Section 2, the sum as defined by Eq. ( p2.37| ) does not exactly 
coincide with the principal value of the Borel integral. The difference is tiny for the 
value of masses considered here and does not affect our conclusions which are based on 
the behaviour of perturbation theory in low and intermediate orders. Second, for the 
bottom quark the BLM scale is defined without taking into account the charm threshold 
in the running coupling. 

The difference between the pole mass and the MS-renormalized mass (normalized 
at mjyfg) for bottom quarks can now be written as 



^ = ^ W+2 ] ( _ # ) asW) 



7T 



Q.91a s (m b ) + higher orders 



(4.17) 



where the second term in square brackets accounts for the genuine gluonic two-loop 
corrections, cf. Eq. ( |4.1| ). Unknown higher order corrections include genuine three- and 
higher loop corrections, vacuum polarization insertions into two-loop corrections as well 
as effects of two-loop running on lowest order radiative corrections. Numerically, for 
ot s {mb) = 0.2, we find 



5m b 
m b 



(16.3 ±2.9 ±1.5)% 



(4.18) 



Compared to the two-loop expression, the estimate obtained from NNA increases the 
mass shift from about 12% to 16%. The second error reflects our estimate of unknown 
higher order corrections, which we allow to be as large as the genuine two-loop correc- 
tions. The first error, which dominates the total uncertainty, represents an estimate of 
the ultimate accuracy of 5m b /m b due to the divergence of the perturbative series. It 
can not be reduced by calculating higher orders (but of course refined - the numerical 
value we quote has been obtained by increasing the uncertainty of M^ +2 \a s ) in Table f| 
by 50%, upon which it is close to the minimal term of the series at n ~ 2 — 3). In 
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absolute values, the intrinsic uncertainty of Sm^ and therefore of the pole mass of the 
bottom quark ranges between 100 and 150 MeV. The above errors do not include an 
uncertainty in a s {m^), which has been fixed to 0.2. 

All results were given in terms of ^ms(^ms)- To run to a different renormalization 
point fi, one should use the anomalous dimension to the same large-/?o approximation. 
In the MS-scheme it is given by 



7m(a s ) 



fi 2 dm(/j,) a s 
m djj? 3tt 



G (-p a s 



(4.19) 



with Gq as in Eq. (fOj) in agreement with [39]. Thus, in this approximation 



m(fi) = m(m) exp 



da* 



37r/3oa, 



Go(-(3 a s 
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J 



(4.20) 



To the approximation considered, one may replace the exponential by the first two 
terms of its expansion. In practice, it might be better to consider the approximation as 
an approximation to the exponent and to keep the exponential. 



5 Scale-setting at next-to-leading order 

So far we have been discussing a scale-setting procedure that utilizes only one-loop 
evolution of the strong coupling. At the same number of loops, where this procedure 
extends the familiar scheme of Brodsky, Lepage and Mackenzie, one also encounters 
diagrams associated with two-loop evolution, which are not suppressed compared to 
insertions of single fermion loops by any small parameter. This Section is devoted to 
the possibility of and difficulty in extending the scale-setting to next-to-leading order 
(NLO). The precise meaning of NLO in this context is illustrated in Fig. [5]: For any 
process, we will calculate the class of higher order corrections, generated by substituting 
the gluon (photon) propagator in lowest order by the chains of Fig. |5]a (this was done in 
previous Sections) and Fig. [5]b. In the abelian theory, the results are exact, while in the 
nonabelian theory, one must again define how to restore the coefficients /5o and (5\ of the 
nonabelian /3-function. Let us note that as in the case of one-loop running this extension 
is heuristically motivated by the behaviour of perturbation theory in large orders. Two- 
loop running is known [|4(], [11]] to modify the strength of renormalon singularities from 
poles to branch points involving the ratio (3\/0q. The class of diagrams in Fig. [5|b is 
expected to dominate at large n over one- loop running by a factor /3i//3q In ri and is 
the first in a series of multiple two-loop insertions that exponentiate to produce an 
enhancement factor n - ^ 1 ^. We are thus again led by the expectation that a class of 
systematically large corrections can be identified and calculated. Consider a physical 
quantity, with perturbative expansion written as 

R-Rtree = 5>„,a s (Q)" +1 , r n = r n0 + r nl N f + . . . r nn N] , (5.1) 

n=0 
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n jwv vwni)vW (no) jwdn^ 



(a) 



(b) 



Figure 5: QED-like diagrams incorporating evolution of the coupling to leading (a) and 
next-to- leading order (b). A circle with letter m denotes a chain of m fermion loops. 
At order a k+1 , the relevant diagrams are specified by n = k and n\ + n 2 + n 3 = k — 2. 
At NLO, the diagrams, where the n2-chain forms a self-energy-type insertion, are not 
depicted. 



assuming for simplicity that R depends only on a single scale Q, which is equal to the 
renormalization scale, such that the r n are numerical coefficients. 

Since the scale-setting prescription is derived from QED (that is, its fictitious version 
with Nf massless fermions), it is useful to consider first the abelian analogue of R. Due 
to the Ward identity, the evolution of the coupling in the abelian theory is generated by 
radiative corrections to the photon propagator. It is then natural to define the scale of 
the coupling at each order in perturbation theory by replacing each photon propagator 
1/k 2 by the full propagator 1/(A; 2 (1 + U(k 2 )) and absorbing the effect of integrating 
loops with full propagators into the normalization of the coupling constant to the order 
where the corresponding skeleton diagram appears. In this way, one is led to a modified 
expansion U 



Specifically, the scale of the coupling at order a s is determined by the replacement 



inside the loop-integration over k, where k is the momentum of the virtual photon. We 
have attached a subscript R to quantities that depend on the choice of the renormaliza- 
tion scheme R. On the right hand side we notice that oty{Q) is the effective coupling, 
defined by the potential between two static sources in momentum space, 



R - Rtree = r a s (Q*) + <WQ**) 2 • • • 



(5.2) 



an(Q) -> a v (k) 



oir(Q) 



(5.3) 



l + U R (k 2 /Q 2 ,a R (Q)) 



V(Q) 



Ait ay (Q) 

Q 2 



(5.4) 
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which defines the V-scheme M. Thus, in any scheme R, the scale Q* R is given by 
averaging the lowest order radiative corrections to the quantity R with the running 
coupling in the V-scheme, ay(k), and notQ with or(/c). By construction an(Q* R ) is 
scheme- independent. The choice of ay(k) to average the loop momentum distribution is 
physically appropriate, because the photon exchanged between static sources at distance 
r has momentum k ~ 1/r. Therefore it is in this scheme that a(k) can be interpreted 
as the effective coupling of a virtual photon of momentum k. 

With this remark one is prepared to adopt the same definition in the nonabelian 
theory. Q* is defined as || 

r a R (Q*) = J d 4 k F(k, Q) , (5.5) 

where F(k, Q) is the integrand of the lowest order radiative correction and 
Cp4:7rav(k)/k 2 the potential between two static colour sources in momentum space, 
given by a Wilson loop. The Wilson loop is nonperturbatively well-defined (up to a 
distance-independent contribution) and the above integral could be evaluated without 
any ambiguity. In the following, we will use it only perturbatively, in the sense of an 
expansion of ay(k) in ay(Q). This restriction is not only self-imposed. Power cor- 
rections to the static potential might be different and in particular larger than those 
to the quantity R, in which case one would not like to use Eq. ( |5.5| ) with power-like 
accuracy. Let us also note that one-loop gauge boson exchange exponentiates exactly in 
the abelian Wilson loop but not in the nonabelian case. Therefore re-expressing ay(Q) 
in terms of some other q:r(Q) inevitably involves contributions to the Wilson loop with 
more than one gluon exchanged. 

The exact (perturbative) evaluation of Eq. ( |5.5| ) is impossible even in the abelian 
theory. One has to resort to some truncation of (a) the /3-function in the ^/-scheme, 
needed to relate ay(k) to a>y{Q) and (b) the relation between the coupling in the 
V-scheme and some other scheme R, if one wishes to express the result in terms of 
aR. To leading order, the /9-function has been replaced by j3 a 2 . After expansion of 
a{k) in a(Q), this corresponds to the diagrams of Fig. ^|a in QED. In this section we 
discuss a truncation, which is guided by the subleading ^"/-dependence of coefficients 
and incorporates the set of diagrams depicted in Fig. ||b in addition to those in Fig. [|a. 
Let us first consider rii = n 3 = in Fig. [|b. The diagram with n 2 = gives flice 3 
to the abelian /3-function. For n 2 > 0, these diagrams give contributions to higher 
order coefficients to the /3-function. Again in QED, this contribution can be written as 
6 n2+ i(— f3o) n2 with a iV^-independent number b n2 . When expanding ay{k) as 

18 At the level of one-loop running, the relevant transformation of schemes is simply a shift of scale, 
since, to this approximation, 
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ay(k) = exp 



a 



a'=a v (Q) ' 



(5.6) 



the diagrams of Fig. |5] are recovered by using 

oo 



(5.7) 



n=l 



but keeping only terms with at most one power of /3i. A different truncation could be 
imagined, where the /3-function is kept only up to fiia 3 but all powers of fli are taken 
into account in the expansion of Eq. (|5.6|) . The difference to the previous truncation 
appears first at order a 5 . 

The transition to the nonabelian theory is arguably uniquely achieved at leading 
order by replacing the abelian (3 by its nonabelian value. At NLO one has to make 
a choice which is less compulsory: Calculate the insertions of the diagrams of Fig. |5|b 
and express them as c n (— /3i)(— Po) n ~ 2 a(Q) n+l in QED. Then replace (3$ and (3\ by their 
nonabelian values. In effect this implies substitution of the true /3-function in the V- 
scheme by Eq. (|5.7|). We can now write the perturbative coefficients of the quantity R 



as 



r n o + 



Pl(-P Q ) n - 2 Cn + (-PoTd n 



n>2, (5. 



where, according to the prescription, the last two terms are absorbed into Q*. Both 
sets of coefficients c n and d n are computed from abelian diagrams. Note that d n is 
constructed as to eliminate the highest power of Nf. On the other hand, to NLO, 
the remainder 5 n still contains subleading flavour dependence, NJ~ , which arises from 
insertion of fermion loops into diagrams with two gluon lines as well as the effective 
three-gluon coupling generated by attaching three gluons to a fermion loop. 

Before turning to the calculation of c n , we want to illustrate the prescription at 
order a 3 . Up to this order, the expansion of R in a certain scheme with coupling a s 
can be written as 



R - R tree = r [a s (Q) + {6 1 + (-/%)<*!> » S (Q) 2 + {<5 2 - P1C1 + (-(3 ) 2 d 2 } a s (Q) 3 ] . (5. 



From Eq. ( |5.5| ) we obtain 

a s (Q*) = a v (Q) + (-/3 )d iav (Q) 2 + + {-Po) 2 d 2 } a v (Q) 3 (5.10) 

where the (scheme-dependent) d n are computed as in Section 2 from diagrams with 
fermion loop insertions. By comparison we find that in the V-scheme {ot s (Q) = otv{Q)) 
c\ = d\. The remaining linear ^/-dependence in 5 2 can then be absorbed into a scale Q** 
of the a^-correction M. Notice, however, that while Q* is defined without knowledge 
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of the exact a^-coefficient, the definition of Q** in the nonabelian theory requires the 
exact af-result. 

To obtain the scale Q* as an expansion in the coupling different from the V-scheme, 
one has to relate the couplings to the same accuracy and expand ay(Q) in, say, «Mg(<5) 
in the same form as Eq. ( |5.9| ): 

aviQ) = «ms(Q) + {71 + (-A)) (-l)} a Ms(Q) 2 

+ {72 -&(-§ + 4C(3)) + (-/?o) 2 y } a m (Q) 3 (5.11) 

One can then determine c% in the MS scheme. Once di, d 2 and c\ are determined, Q 2 
to a^-accuracy is given by 

Ql = Qe- d ^(^ + {{-fo)b 1 + b 1 ^a 8 {Q)} . (5.12) 
bi and 61 are uniquely determined by the condition 

a s (Q* 2 ) = a s (Q) + (-/? )<WQ) 2 + {(-/3o) 2 (d? - 260 - (di - 2&i)fr} a s (Q) 3 + . . . 

= a s (Q) + (-/3 )dia(Q) 2 + {(-/5o) 2 d 2 -/?iCi}a s (Q) 3 , (5.13) 



which specifies that the diagrams of Fig. |5] relevant at order a 3 are absorbed into Q 2 in 
any scheme. In the V-scheme b\ = 0, but in general b\ is non-zero. 

In the V-scheme the prescription outlined here coincides exactly with the one by 
Brodsky and Lu [0. When Q\ is defined in other schemes, we suspect that Q* defined 
in does not absorb the QED diagrams of Fig. |5|b and part of the flavour-dependence 
from these diagrams is hidden in Q**, since 61 is always zero in [0. 

When we compute c n in the following Section, rather than using the V-scheme in 
intermediate steps, we will perform the average of Eq. ( |5.5| ) directly in an arbitrary 
scheme by replacing av{k) with the appropriately truncated form of Eq. ( |5.3|) . 



5.1 Abelian diagrams at NLO 

We now evaluate the perturbative coefficients generated by inserting the chains of Fig. [5] 
into the gluon (photon) line which appears in lowest order radiative corrections to an 
observable R. In this subsection (3q and (3\ refer to their abelian values: (3q = Nf/ (3ir) 
and fix = Nf/(4iT 2 ). These chains have previously been considered for particular quanti- 
ties: The muon anomalous magnetic moment |16| and the self-energy of a static quark 
ffT| . Because vacuum polarization insertions are universal, the results can easily be 
adapted to the general case. The vacuum polarization to the present approximation is 
given by the diagrams with n = 1 and n\, n% = in Fig. |5| (Q 2 is euclidian): 



-Aj«(m) 




(5.14) 
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When acv(k) is expanded inside the integrand of Eq. (|5.5|) , the result is expressed 
in terms of coefficients of the vacuum polarization and averages of ln n (k 2 /Q 2 ) in the 
lowest order radiative correction F(k,Q), which have already been evaluated through 
the insertion of the diagrams of Fig. [5]a. Therefore, since 111 is known, the extension to 
incorporate Fig. |^b is merely a combinatorical problem. To organize the combinatorics, 
it is convenient to introduce the Borel transform. First, we define the leading order 
and next-to-leading order Borel transform as the result of insertion of diagrams with a 
single chain (Fig. |5|a) and the sum of both diagrams in the Fig. |5|. From Eq. ( |2.6| ) and 
Eq. (|5j): 



OO j 

B LO [R](u) = r £-V 

n=0 
oo 



jNLO a 

B NLO [R](u) = ro^^-n" d ^° = d n -^c n (n > 2) (5.15) 

n=0 U - 



and d^ LO = d n for n = 0, 1. After Borel transformation of Eq. ( |5.5| ), the truncation of 
av(fc) discussed above corresponds to insertion of 



B 



a(Q) 



l + n(A: 2 /Q 2 ,a(Q)) 



(u) 



Q 2 



Pi 
^6 



B 



Hi 



(u-v) (5.16) 



into Eq. (|5.5| ) instead of the complete ay(fc). On the right hand side, we have neglected 
multiple insertions of III. The Borel transform of Hi /a can be represented as 



B 



Hi 



a 



F(u) - G(u) 



(5.17) 



where F(u) is a scheme-independent function that can be obtained from Eq. (|3.9| ) by 
integration with respect to Q 2 : 



F(u) 



32 



3 l-(l--u) 2 



E 

k=2 



{-ifk 



(k 2 



l-u 



,2\2 



(5.18) 



n=-l 



G(u) is a scheme- dependent integration constant. One can use the renormalization 
group equation obeyed by the photon vacuum polarization to relate the expansion 
coefficients of G(u) to those of the /3-function in the chosen scheme. The precise relation 
is as follows: Let us write the highest power of f3 n as 



n\m 



KPi(-Po 



,n-l 



(5.19) 



Then 
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g{u) = y: 



b n 



+l u n-l 



9n 



K 



+2 



(n + 1)! 

With these preliminaries, we insert Eq. (|5.16 ) into Eq. ( |5.5|) and obtain 



(5.20) 



B Nh0 [R](u) = B LO [R](u)-% [ — (B LO lR}(u)-B L o[R}(v)) 

J u-v 



u 

~%\ dvv(B LO [R}(u)F TCg (u -v)- B LO [R](v)G Teg (u - v)) , 
p J 



(5.21) 



where F reg (u) and G reg (u) are defined by removing the pole at u 
F{u) - l/u, G rcg (u) = G{u) - 1/u. 
The coefficients <i^ LO are given by 



0: F reg (u) 



1 d 



^NLO 
n r du 



•it 



-Bnlo[R](u)\ u = ■ 



(5.22) 



Taking the derivatives of Eq. (|5.21|) we obtain (for n > 2) 



-NLO 



0i 

d ™ ~ IM 71 (^( n + 1) - V>( 2 )) d n -i 
Po 

n-2 

fn-2-k ~ (k + l)g n -2-k 



p o fc=0 



n 
k 



(n-2-k)\d k , (5.23) 



where ip(x) is the logarithmic derivative of the T-function. Note that for large n the 
second term on the right side indeed dominates the first one by a factor p\//3 2 In n. With 
correspondingly changed conventions, this equation agrees with the corresponding one 
in [16| . The expansion coefficients of F(u) are given by (16 



fn 



2, 

-g(n + 2) 



n + 7 16 
-2n - 2 — H 

2 n+3 ^ n + 2 



^ fc(l-2- 2fc )(l-2 2fc - n - 3 )C(2A: + l) 



fc=i 



(5.24) 



where £(fc) = Y^?=i n k an d [••] denotes the integer part of the number in brackets. 
The coefficients g n depend on the scheme employed for the definition of a. In the MS 



scheme, we can use the /3-function to the approximation required here from [ 39| , [T6 
We then find 



9™ 



1 



d n+l 



(n + l)!(n + 2)! du n + l 



[1 - u)(l + 2u)(3 + 2tt)r(4 + 2u) 

9r(2 + u) 2 r(3 + u)r(i -u) 



(5.25) 



\u=0 
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In the V^-scheme, we have 



9n — fn ■ 



(5.26) 



The simplest way to see this is that in this scheme by definition 11(1, ay{Q)) = 0. Then 
the expansion of the integrand in Eq. ( |5.5| ) has only non-zero powers of ln(k 2 /Q 2 ) and 
d can not appear in Eq. ( |5.23|) . 

5.2 Numerical analysis 

The transition to QCD according to the prescription formulated above is performed by 
replacing (3 2 in Eq. (|5.23|) by its QCD value. Note that contrary to d n the NLO 
coefficient depends explicitly on the number of flavours. In this Subsection, rather 
then presenting values for the BLM scale Q* to NLO accuracy, we give the perturbative 
coefficients in low orders generated by expansion of of a s (Q*) in the MS-coupling, that is 
the coefficients d^ LO and M^ LO (— Poa s ), defined in complete analogy with M n (— /3oa s ). 
To compare the coefficients obtained from keeping only vacuum polarization with the 
exact ones, if available, we shall also define c/^ xact as the exact coefficient, divided by 
( — Numerical values for the difference between pole and MS mass are shown in 
Table |5|, and for the derivative of the hadronic vacuum polarization Q 2 dU/dQ 2 and the 
r decay width in Tables ^| and 0. 

The main conclusion for the radiative mass shift is that the effect of two-loop running 
remains small up to the order, at which the series has to be truncated. It leads to a less 
then 5% additional increase of the one-loop radiative correction, which can be neglected 
in view of the uncertainties that have been discussed in Sect. 4. This is reassuring, 
because it supports the suggestion that the mass shift is to a large extent given by one- 
loop running effects, although it must be borne in mind that by incorporating two-loop 
running one does not gain control over genuine two- and higher loop corrections. 

On the other hand, the effect of two-loop running is very large for the derivative of 
the vacuum polarization and the efficiency of extending the BLM prescription beyond 
leading order can be doubted in this case. Although the inclusion of two- loop running 
improves the estimate of g?2 as compared to the exact coefficient d^ (Table |^), for n > 3 
the effect of two-loop running is exceedingly large and this improvement might as well 
be accidental. 

It is interesting to observe that for the r hadronic width itself (Table |7|), such irreg- 
ularities do not occur, which is in qualitative agreement with the conclusion of Sect. 3.1 
that the perturbative expansion of R T has a smoother behaviour in intermediate orders 
of perturbation theory than Q 2 dU/dQ 2 . The relative importance of two- loop running 
increases with n as expected from the asymptotic Inn-enhancement. Though smooth, 
the effect of two-loop running is significant for R T and points towards an even further in- 
crease of the cumulative effect of higher order perturbative corrections and consequently 
a further decrease of a s (m T ) from R T . 
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n 




e° 


M n [b-quark] 


M^ LO [b-quark] 





1 


1 


1 


1 


1 


4.68615 


4.68615 


1.623 


1.623 


2 


17.6227 


19.6884 


1.935 


1.972 


3 


109.859 


127.529 


2.193 


2.272 


4 


873.924 


1138.28 


2.467 


2.628 


5 


8839.69 


12085.6 


2.835 


3.131 



Table 5: Comparison of leading order and next-to-leading order coefficients for the 
difference of pole and MS mass, where quark masses inside loops are neglected and 
loops containing the heavy quark are excluded (cf. Sect. 4). The comparison of the 
modification of lowest order radiative corrections is for a s (rrib) = 0.2, relevant to bottom 
quarks. 



n 




c° 


jexact 


M n 


M NLO 
n 


?M~exact 
n 





1 


1 


1 


1 


1 


1 


1 


0.691772 


0.691772 


0.728809 


1.158 


1.158 


1.167 


2 


3.10345 


2.31417 


1.25847 


1.321 


1.271 


1.233 


3 


2.18004 


7.40378 




1.347 


1.360 




4 


30.7398 


18.4580 




1.432 


1.411 




5 


-34.5336 


146.293 




1.410 


1.503 





Table 6: Comparison of leading order, next-to-leading order and exact coefficients 
of Q 2 dU(Q)/dQ 2 . The NLO and exact values are given for Nf = 3 and the partial 
sums M n for a s = 0.32. Exact values, when available, include genuine higher order 
corrections. 



n 


dn 


jNLO 
a n 


jexact 


M n 


M NLO 


?M~exact 
lvl n 





1 


1 


1 


1 


1 


1 


1 


2.27511 


2.27511 


2.31213 


1.521 


1.521 


1.529 


2 


5.68475 


5.98780 


5.20806 


1.819 


1.835 


1.803 


3 


13.7536 


18.1248 




1.984 


2.053 




4 


35.1470 


54.7939 




2.081 


2.203 




5 


84.4066 


178.897 




2.134 


2.316 





Table 7: Same as in the previous table for the perturbative expansion of the r hadronic 
width (a s (m T ) = 0.32). 
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In Appendix C, we give exact results for the five-loop diagrams that enter the 
calculation of d^ LO for Q 2 dU/dQ 2 . 



6 Conclusions 

In this paper we have shown how to deal with higher order vacuum polarization in- 
sertions into radiative corrections to observables in QCD. This enterprise has been 
motivated by the fact that these higher order corrections often lead to disturbingly 
large coefficients in the perturbative expansion. This trend is systematic, but since it 
can be identified, we suggest to calculate these corrections and separate them from the 
remaining ones. The behaviour of the remaining series should then be improved. Even 
if the remaining corrections turn out not to be small, the summation of vacuum polar- 
ization insertions can be motivated as a physical way of scale-setting. In this respect 
it can be considered as an extension of the scale-setting proposed by Brodsky, Lepage 
and Mackenzie |TJ and Brodsky and Lu [0, as far as two-loop running is concerned. It 
is worth noting that the familiar BLM scale typically overestimates the size of higher 
order vacuum polarization insertions. The smallness of this scale needs not necessarily 
indicate a failure of perturbation theory, unless it is indicative of its divergence already 
at two-loop order. 

We would like to stress the computational ease with which the summation of effects 
of one-loop running can be performed compared to a genuine higher order calcula- 
tion. We thus have devoted considerable space to technical aspects, which should allow 
routine implementation of this summation. The two existing techniques can be summa- 
rized as follows: The first one [T3, IIR O, |3[ amounts to direct evaluation of the relevant 



higher order Feynman diagrams. A convenient tool to organize such a calculation is 
the Borel transform. The result can be obtained as a certain analytically regularized 
Feynman integral plus - if necessary - a subtraction function that is easily computed 
directly from a large-mass expansion. This technique has its limits, when applied to 
processes with several scales or to physical cross sections that are more directly cal- 
culated as a sum of real and virtual corrections. In these cases, we have applied a 
dispersion technique [L^, |J (or, if one starts from the Borel transform, a Mellin 
transformation) to reduce the problem to calculation of lowest order corrections with 
finite gluon mass. It is in this representation, that summation (up to an accuracy set 
by renormalons) is most easily performed. This summation is usually complicated by 
the need to analytically continue the Borel transform beyond its radius of convergence 
and to take the (principal value) integral over the Borel parameter. Our Eqs. ( |2.24| ) 
and ( 2.29|) are immediately suited to numerical evaluation. Given this simplicity, the 



computational expense seems worth the investment even if one could only hope to ab- 
sorb a higher order radiative correction that has the correct sign compared to the exact 
one. 

In the present paper, we have investigated two quantities in more detail. Higher or- 
der /3q a™ +1 -corrections are indeed sizeable for the hadronic decay width of the r lepton, 
when expressed as a series in a s . We compared fixed order perturbative approximations 
to R T with those from a partial resummation of running coupling effects due to contour 
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integration and found that in intermediate orders the first one has a smoother be- 



haviour, since potentially large corrections from contour integration conspire with the 
divergent coefficients of Q 2 dYi/dQ 2 to produce an effective suppression of ultraviolet 
renormalons (Table [l]). Resummation of all one-loop running coupling effects leads to 
an approximate 10% decrease of a s (m T ) determined from hadronic r decays to a central 
value o!0 



a s (m T ) ~ 0.29. 



(6.1) 



This shift of about one error margin of previous analyses is not caused by a single 
large next-order coefficient, but by the constructive addition of several higher order 
coefficients beyond the exactly known coefficient at order a? s . 

The accurate extraction of a s from r-decays relies on the absence of AQ CD /m^- 
corrections. Although from a purely theoretical point of view, the situation with respect 
to such terms is not conclusive, given that the limitations of duality are only poorly 
understood, we conclude that effects associated with summation of the divergent per- 
turbative expansion should be excluded as a source of such terms, provided one accepts 
their absence in the relevant current correlation functions in euclidian space. 

For the difference between the pole and MS-renormalized mass of a heavy quark, we 
conclude that the large two-loop correction found in is probably not accidental but 
the first one in a rapidly divergent series of higher order corrections dominated indeed 
by one-loop running of the coupling. For the charm quark, the divergence is such that 
the perturbative series can certainly not be improved beyond two-loop order. For the 
bottom pole mass, our estimate incorporating one-loop running reads 

^ = (16.3 ±2.9 ±1.5)%, (6.2) 

which is about 30% larger than the two-loop result without resummation. We expect 
this difference to be significant in phenomenological applications. The first error quoted 
is associated with the divergence of the series that relates the two mass definitons and 
is irreducible. Numerically this error amounts to an uncertainty somewhat larger than 
100 MeV for the mass difference, a value in between the estimates initially reported 
in |37| Short-distance observables containing quark masses, when expressed in 

terms of pole masses will generally exhibit large coefficients implanted by the use of 
this parameter. Thus, although both mass definitions are ultimately unphysical, we 
expect that the MS mass can be determined to better accuracy through perturbative 
relations^]. 

Finally, let us mention that the calculation of vacuum polarization can never replace 
an exact calculation. Until this is completed, we find it worthwhile to incorporate the 
systematic effects exhibited by one-loop running. Surely, the authors are among those 

19 As an input we have taken R T — 3.56 ± 0.03 from |2^| , see Sect. 3.1. 

20 This argument does not apply to pole mass differences. The large coefficients associated with 
the pole mass can be understood as a universal (flavour-independent) additive mass renormalization, 
which cancels in the difference. 
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who await further exact results with suspense. 



Note added. While this paper was written, we were informed by M. Neubert on 
work of his, which partially overlaps with the present one. We acknowledge exchange 
of manuscripts and are very grateful for the ensuing discussions which helped to clarify 
our presentation. We disagree with the conclusion of (|2]| that the difference of various 



summation prescriptions for cross sections should be interpreted as a sign and quan- 
tification of the failure of the operator product expansion in the physical region. As 
emphasized in Sect. 3.2, once one abandones Borel-summation type prescriptions, which 
are distinguished by their relation to the OPE, one can introduce l/Q 2 -differences in 
euclidian space just as well and there is no discrimination between the euclidian and 
minkowskian region from the point of view of summation of perturbation theory, given 
the present state of knowledge. 

Acknowledgements. M. B. would like to thank D. Broadhurst and A. Pich for 
stimulating discussions. M. B. is supported by the Alexander von Humboldt-foundation. 
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Appendices 



A Subtractions 

In general, we will also be interested in ultraviolet divergent quantities. Examples 
include the difference between the pole mass and the MS-renormalized quark mass 
and the correlation function of heavy-light currents in heavy quark effective theory. 
Denoting a generic quantity by R(a), the renormalized Borel transform has the form 

B[R)(u) = B[R] (u) + S R (u) , (A.l) 



where B[R] (u) is the bare Borel transform, calculated from diagrams without overall 
subtraction, using the Borel parameter itself as a regulator. This corresponds to ana- 
lytic regularization |l5j . Consequently, the bare Borel transform is singular at the origin 



(Recall that the perturbative series is generated by derivatives of the Borel transform 
at the origin.). This singularity is subtracted by the function Sr(u), which depends 
on the renormalization scheme. In minimal subtraction schemes uSr(u) is an entire 
function (at least for one chain of fermion loops) and - when combined with the ap- 
propriate singular piece of B[R] (u) to be defined below, Eq. ( |A.17| ) - Sr(u) yields an 
unambiguous contribution to the bubble sum 

oo 

Sfl(a) = ("^-) J due- u ^-^ (Sr(u) + singular term in B[R] Q {u)) . (A.2) 
o 

General expressions for Sr(u) in minimal subtraction schemes have been derived in 
Appendix A of JT7J. In this Appendix we outline a simple method to obtain Sr(u) 
and calculate Sr(cx) in minimal subtraction schemes. An essential simplification comes 
from the fact that one needs only the dominant large-A behaviour of dimensionally 
regularized Feynman amplitudes, where A is a mass for the gluon. As an illustration of 
the technique we calculate the subtractions for the correlation function of two heavy- 
light currents in heavy quark effective theory. 



A.l Calculation of Sji(a) 

As always we assume that the quantity R has been made dimensionless, is infrared 
finite and such that the order-a radiative correction comes from gluon exchange. We 
shall also assume that the corresponding diagrams have one loop and are at most 
logarithmically ultraviolet divergent. The last assumption is actually unnecessary, see 
the explicit example in the subsequent subsection. The dimensionally regularized order- 
ed correction can then be represented as 

r^ e (e) a = f d d k F{k, Q)-L . (A.3) 
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We use d = 4 — 2e and denote by Q a set of external momenta. For gauge-dependent 
quantities, we take Landau gauge. The factor (k^k u /k 2 —g^ u ) from the gluon propagator 
is included in the integrand F(k, Q) and we do not write Lorentz indices explicitly. We 
define the both, dimensionally and analytically regularized coefficient by 

f bare ( S ,e)a = a^J d d kF(k,Q)± (~^\ , (A.4) 

where \i is a subtraction scale. Note that B[R]q(u) = fQ are (-u, e = 0)e~ uC , see Eq. ( g3g ). 

Insert n fermion loops into the gluon line of all diagrams that generate r . The 
fermion loop integrations can be done. Each loop gives a factor 



/3 / 6r(l + e)r(2-e 



,2 



e T(4 - 2e) \ 4vr^ 2 y ' 

where (3q = T/(3tt) is the fermionic contribution to the beta-function (T = 1 in QED 
and T = 1/2 in QCD). Performing the final integration over gluon momentum k, the 
result for the coefficient of order a n+1 can be written as 



„bare 



r- 



(e)= (n + l)(-e)n + l G(-e,-(n + l)e) (A.5) 



G(-e, -(n + 1)6) = ^(4^)- 6r(1 + 4 e) _ r( 2 2 g) e) j (n + l)(- e ) f bare (* = nc, e) , (A.6) 

where Eq. ( |A.4j ) has been used. The function G introduced above coincides with the 
one of Appendix A of |T7| [Note that d = 4 + 2e has been used in |T7|], which motivates 



the signs in the arguments of G above.]. To obtain the subtraction function Sr(u) in 
the MS scheme, one has to add to r bare (e) the diagrams with fermion loops replaced by 
their MS counterterms and then compute the finite part of the sum. For any quantity 
R, the steps of Appendix A in [T7J for the self-energy of a heavy quark can be repeated 
with the general result 

S R ( U ) = ^1 G (u) = f;^u», (A.7) 

n=0 

where g n are the expansion coefficients of 

oo 

Go(-e) = G(-e, -(n + l)e = 0) = £ g m {-e) m . (A.8) 

m=0 



Since Go is related to subtractions, it should be unnecessary to calculate the full 
diagram, encoded in the function G, in order to obtain Go. We prove that Go can be 
deduced from the large-mass expansion of the dimensionally regularized one-loop coef- 
ficient rQ are (A 2 , e) computed with a gluon mass A. That is, 1/k 2 in Eq. ( |A.3| ) is replaced 
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by l/{k 2 — A 2 ) (in Landau or Feynman gauge). Since the original four- dimensional 
integral was logarithmically ultraviolet divergent, the asymptotic behaviour at large A 
is given by 

r w (A2)e) A^oo_i roo(e) ^y _ (A9) 



e 



This equation defines the A 2 -independent function r oc (e). Its value at e = is denoted 
by r^. By definition 

G »<- £) = -jhymr^w (s + £) f ° b "' (s ' £) ■ (A ' 10> 

To evaluate the limit, we use that fo are (s, e) is related to the one-loop coefficient with 
finite gluon mass rQ are (A 2 , e) by a Mellin transform: 

rr e (s, e) = J ^ (4) S -o arC (A 2 , e) (A.ll) 

^ o ^ ' 

This represents the analytic continuation of the analogon of Eq. ( |2.17| ) into the u- interval 
] — e, 0[. Integrating by parts, we get 



x "\ ~(*+e) 

vr (s + e) J v /i 2 



r^%s,e) = ld\ 2 (-A f'{\\e) , (A.12) 



where /(A 2 , e) is defined by dividing fo are (s, e) by (/i 2 /A 2 ) e . The limit is now easily 
taken with the result 

lim (s + e) f bare ( S , e) = roo (e) . (A.13) 

This yields the final result 

_ , , 1 r(4 + 2u) sin7rn , . , . „ 

(47r) " 61 (1 — ujl (2 + u) z iru 

The computation of the subtraction function has been reduced to calculation of the 
large mass limit of one-loop corrections with non-zero gluon mass. 

Combining Eq. ( |2.17| ) for the bare Borel transform with Eq. ( |A.7| ), the renormalized 
Borel transform is given by 

. sin7rn f „, / A 2 ^\ Gn(u) 
B[R)(u) = / dA 2 e c r (X ) H . A.15 

TTU J \ LL I U 

V 7 

The derivative ^(A 2 ) is ultraviolet finite, which allows to put e = 0. However, the 
integral is not yet finite at it = and the above expression is not suited to take 
derivatives at u — 0. Since, by Eq. ( |A.9| ), Tq(A 2 ) = r^/A 2 for large A 2 , we obtain 
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B[R](u) 



smnu 



TTU 



d\ 2 



A 2 



e(A 2 



+- G (u) 
u 



Sill TTU 



TTU 



(A.16) 



In this form the integral exists for u = and the second line is finite at u = 0. We note 
that this expression is equivalent to, but slightly different from the one given in |4j]. In 
the present form /i 2 and C appear only in their natural combination n 2 e~ c . 

To obtain the bubble sum, all steps that lead from Eq. fl2.17| ) to Eq. ( |2.24| ) can now 
be repeated on the first line of Eq. (|A.16|) . The Borel integral of the second line is given 
by 



oo 

S R (a) = (-i-) J due~ u /^ a ^(G (u)-r c 
o 



sm7rn 



TTU 



-uC 



(A.17) 



-foot 




arctan(7r/3 a) 1 / 2 2 2 \ 

a + -In 1 1 +vr /3 a 1-1 

7rpo« 2 V / 



Thus, in case subtractions are required beyond coupling renormalization, Eq. ( |2.24| ) is 
replaced by 



r a s Moo(a s 



+ 



d\ 2 d>(A 2 ) (r^(A 2 ) - ^ 6(A 2 - ^ 2 e" c )) + [r (A 2 ) 



ro(0)] 



du 



{G (U) - Too) + 

' c 



arctan(7ra s ) 1 , / 9 9 , 

+ - In (1 + vr 2 a? ) - 1 

7ra, 2 V 



(A.18) 



with G (u) given by Eq. (|A.14|) . 



A. 2 A sample calculation 

As a non-trivial example, we calculate the subtraction function for the correlation 
function of heavy-light currents in heavy quark effective theory, 



n 5 M = i / d 4 *^*) (0\T{jl(x)j 5 (0)}\0) j 5 (x) = h v (x)i 75 q(x) . 



(A.19) 



The corresponding bare Borel transform has been given in [17|. At first sight the 



method exposed in the previous subsection appears inapplicable, because the diagrams 
to be considered have two loops, see Fig. ^, and, because the correlation function is 
quadratically divergent, one also needs subtractions for the two-loop diagram itself. To 
eliminate these, we shall consider the third derivative 
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a 




Figure 6: Radiative corrections to the correlation function of heavy-light currents. 
The shaded circle represents a current insertion with momentum q and the double line 
denotes a heavy quark propagator. 



D{ U ) = . . (A.20) 

Then all subtractions originate from divergent one-loop subdiagrams, which are subse- 
quently inserted in the lowest order one-loop diagram for H^u). We remind the reader 
that the heavy-light current in the effective theory is not conserved. We treat the three 
diagrams in turn and take Landau gauge. 

Let us start with diagram (a) and investigate the one-loop correction to the heavy- 
light vertex with non-zero gluon mass. The numerator of the corresponding Feynman 
integrand is proportional to 

k) lh v T {k 2 9frr - k p k T ) , (A.21) 

where v is the velocity of the heavy quark and the momentum assignments to the lines 
of the diagram are evident. Since the contribution proportional to A _2e required by 
Eq. ( A. 9 ) comes from k ~ A ^> q, we can drop in the above expression and the 



integral is simplified to 



d d k i> k — v ■ k . , . 

(A.22) 



(2Tr) d (q - k) 2 v {q r - k) (k 2 - A 2 ) 
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The corresponding integral with numerator replaced by k^ has structure FiV^ + F 2 q fl + 
F 3 q'^. By dimensional counting, only F 1 can behave as X~ 2t for large A 2 , the other struc- 
tures being suppressed by powers of A. But then the potential term \~ 2e is proportional 
to $ i> — v 2 = and therefore r 00(e) vanishes identically for this diagram (This property 
is specific to Landau gauge.). 

For diagram (b) we need the self-energy insertion of a light (massless) quark, given 

by 



d d k 



(2Tr) d (p - k) 2 (k 2 - A 2 ) 



9pr 



(A.23) 



where p is the external momentum for the subdiagram. Again the term involving A in 
the large A 2 -expansion comes from k ~ A p, which allows to expand the propagator 
— k) 2 . The integrand simplifies to 



(p - k) 2 (k 2 - A 2 ) 



9pr 



kpkj- 



k 2 



e (3 - 2e) 



2-e k 2 (k 2 -X 2 ) 



+ 



(A.24) 



where the neglected terms do not produce a A 2e contribution. The result is 

'47T/x 2 V C F a (3-2e)r(l + e) 



A 2 



4vr (l- e )(2-e) 



(A.25) 



Integration of the final quark loop is straightforward, and taking three subtractions we 
obtain 



G oH u ) 



C F N C 

47T 3 



6 l J r(i -u)r(2 + u ) 2 r(3 + n) 



(A.26) 



Turning to the self-energy of a heavy quark in heavy quark effective theory in di- 
agram (c), we encounter another apparent obstacle to applying the technique of the 
previous subsection: The subdiagram is linearly ultraviolet divergent, and - at large A 2 
- behaves like A 1_2e . The remedy is simply to ignore this term and pick up the term 
A~ 2e , since power divergences are discarded in minimal subtraction schemes. The large 
mass expansion of the heavy quark self-energy is 



H) 



l+jf* / Att/j, 2 ] C F a 



X 2 



47T 



(A.27) 



Discarding the first term and inserting the remaining one into the quark loop, we get 



GS?>(«) 



CfNc 

47T 3 



\ (o 1 2u) F(4 + 2n) 



(A.28) 



Adding the contributions from all diagrams, the result is 
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G (u) = 



CfNc 

47T 3 



2u) 



T(4 + 2u) 
r(l -u)T(2 + u) 2 T(3 + u) 



(A.29) 



B Radiative corections to R e + e and R T with finite 
gluon mass 

We write the lowest order radiative corrections to the total cross section of e + e~ anni- 
hilation and the hadronic r decay width with finite gluon mass A as 



1 + a s (s) {r e T*- (y) + 6(1 - y) r%*- (y)}] , 
1 + a s (m r ) {rf rt (y) + 6(1 - y) r^(y)}] , 



(B.l) 
(B.2) 



where y = A /s or y = A /m 2 ., respectively. Quarks are taken as massless. r (y) are 
virtual and r real (y) real gluon emission corrections: 



„virt 



(v) = 



2 

3vr 



2(1 + y) 2 1 - iln 2 y + lnyln(l + y) + ^+Li 2 (-y) 



- - 2y - 3 In y - 2y In y 



2 

3^ 



2(1 + y) 2 1 X - In 2 y - 2 In y ln(l + y) - ^- - 2Li 2 (-y) 



+ 5 - 5y 2 + 3 In y + 4y In y + 3y 2 In y 



(B.3) 



(B.4) 



and 



rf*(y) = 



- 2577 + 72tt 2 + (2392 - 240vr 2 ) y + (828 - 432vr 2 ) y 2 + 144(1 - vr 2 ) y 3 
-24vr 2 y 4 - 1332 In y + (2208 + 288vr 2 ) y In y + 792y 2 In y + 144y 3 In y - 216 In 2 y 
+720y In 2 y + 1296y 2 In 2 y + 432y 3 In 2 y + 72y 4 In 2 y + 864y In 3 y 



+ (432 - 1440y - 2592y 2 - 864y 3 - 144y 4 ) (lny ln(l + y) + Li 2 (-y)) 

+1728ylnyLi 2 (-y) -3456yLi 3 (-- 

V y 



(B.5) 



r^(y) 



324vr 



2901 - 72tt 2 + (-9736 + 240vr 2 + 5184C(3)) y + (6120 + 432^ 2 ) y 2 



+(456 + 144vr 2 ) y 3 + (259 + 24vr 2 ) y 4 + 1332 lny - (1776 + 864vr 2 ) y lny - 4176y 2 lny 
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(IVl) 



(IV2) 



(VI) 




(V2) 



(V3) 



Figure 7: Representative abelian four- and five-loop diagrams for the hadronic vacuum 
polarization. It is understood that all self-energy type diagrams are added to each class 
of diagrams. 



-288y 3 In y + 216 In 2 y - 720y In 2 y - 1296y 2 In 2 y - 864y 3 In 2 y - 144y 4 In 2 y 
-1440y In 3 y + (-864 + 2880y + 5184y 2 + 1728y 3 + 288y 4 ) (hxy ln(l + y) + Li 2 (-y)) 



C Abelian five-loop diagrams to the hadronic vac- 
uum polarization 

In this Appendix we collect the abelian five-loop diagrams to the hadronic vacuum 
polarization that are included in the scale-setting in next-to-leading order as described 
in Sect. 5. For completeness, we also list the four- loop diagrams. Numerical values are 
given for 



The diagrams are shown in Fig. [7[ Each diagram shown stands for the class of diagrams 
obtained, when one adds the diagrams where gluons form self-energy type insertions 




(B.6) 




(C.l) 
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rather than exchange type topologies. Diagram (V2) also includes the symmetric dia- 
gram, where the two inner loops are interchanged. We obtain: 



7-)( IV1 ) 



D (IV2) 



(V2) 



D 



D (V3) 



C F {TN f f [i§ 
C F {C F TN f ) 

C F {TN f f 
C F {TN f ){C f TN f ) 
C F (TN f )(C f TN f ) 



C(3) 



64" + 2 C(3) 



6131 203 A . o . 5 A/r . 
W + ^C(3) + -C(5) 



3571 59 



576 
10199 



C(3) + 2C(3) 2 



345¥-2^ + C ^ 2 



(C.2) 
(C.3) 

(C.4) 
(C.5) 
(C.6) 



Here T = 1/2, C F = 4/3 for SU(3) with fermions in the fundamental representation 
and T = 1, C F = 1 for U{\). The sum of these terms yields d^ 1 " and c?3 LO in QED 
(after proper normalization). 
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